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was  Program  Manager 


ABSTRACT 


This  document  reports  the  results  of  a  twelve  month  research  study  of  a 
transpiration  cooling  system  concept  having  application  to  the  nose  tip  and 
control  surface  leading  edges  of  an  advanced  ABM.  This  concept  uses  a 
reactive  gas  coolant  and  a  hot  wall  matrix.  Testing  has  been  conducted  to 
determine  the  permeability  and  inertial  resistance  coefficients  of  candi¬ 
date  refractory  porous  and  discrete  hole  matricies.  Testing  was  also  con¬ 
ducted  to  determine  the  amount  of  heat  transfered  from  the  hot  matrix  to 
the  coolant  gas.  The  results  show  that  system  weight  savings  would  be  real¬ 
ized  from  use  of  the  CONAP  concept  when  compared  to  a  water  system.  The 
data  taken  have  been  correlated  with  curve  fits  and  used  to  conduct  prelim¬ 
inary  concept  design  analyses. 
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SUMMARY 


This  investigation  consists  of  the  initial  task  in  the  development  of 
efficient  transpiration  cooling  technology  for  Advanced  Interceptor  Require¬ 
ments  in  order  to  provide  a  lighter  weight,  more  erosion  resistant  system 
than  conventional  water  systems.  The  Controlled  Atmospheric  Protected 
System  (CONAP)  is  a  concept  which  utilizes  a  reactive  gas,  ammonia,  as  the 
coolant  and  a  refractory  metal  as  the  matrix.  The  injection  of  ammonia 
into  the  boundary  layer  medium  cools  the  matrix  and  reacts  with  the  boundary 
layer  oxygen.  In  principle  this  control  of  oxygen  allows  operation  at  high 
surface  temperature,  which  greatly  lowers  the  coolant  flow  requirements 
while  limiting  the  oxidation  of  the  hot  surface  materials  to  acceptable 
values . 

The  purpose  of  this  initial  study  was  to  evaluate  available  porous  re¬ 
fractory  materials  and  their  applicability  as  nose  tip  c.  air  vane  leading 
edge  materials  in  a  flight  environment.  Characteristics  of  the  materials 
and  the  operating  efficiency  of  the  concept  were  evaluated  by  laboratory 
testing. 

These  laboratory  tests  consisted  of  obtaining  flow  characteristic  data 
(permeability  and  inertial  resistance  coefficient)  for  a  variety  of  porous 
tungsten  and  discrete  hole  tungsten  matricies;  also,  the  heat  transferred 
from  a  hot  tungsten  matrix  to  the  ammonia  gas  passing  through  was  measured. 
The  materials  and  concepts  that  were  evaluated  and  the  tests  performed  are 
shown  in  Table  I. 

The  tests  were  performed  on  each  of  the  materials  and  concepts  shown. 
Heat  transfer  coefficient  data  and  flow  characteristics  data  were  obtained 
for  a  wide  range  of  porous  tungsten  materials  ranging  from  50  percent  (by 
volume)  tungsten  to  80  percent  tungsten.  Complete  correlations  of  the  data 
have  been  obtained.  The  results  of  these  tests  and  correlations  are  briefly 
summarized  in  Table  II  and  described  below.  Local  mass  flow  in  porous  tung¬ 
sten  materials  can  be  greatly  increased  to  meet  local  flow  requirements  by 
augmentation  with  blind  holes  without  significantly  affecting  the  internal 
heat  transfer  coefficient.  This  concept  has  potential  usage  for  the  sonic 
point  on  a  nose  tip  or  the  stagnation  line  on  an  air  vane  leading  edge. 
Schlieren  photographs  show  that  the  blind  holes  do  not  starve  other  areas 
of  coolant.  Surface  vapor  deposition  modification  to  porous  tungsten  matri¬ 
cies  can  decrease  the  mass  flow  in  a  local  area  without  affecting  the  heat 
transfer  coefficient.  This  concept  can  be  effectively  used  on  the  sidewall 
of  a  nose  tip  or  air  vane  leading  edge  to  minimize  the  variation  between 
actual  and  ideal  coolant  mass  flow  rates. 


. . 
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TABLE  I 


The  heat  transfer  coefficient  and  flow  characteristic  data,  once  cor¬ 
related,  were  used  to  optimize  the  materials  and  concepts  for  use  in  a 
flight  environment.  The  efficiency  of  the  CONAP  concept  can  be  compared 
to  equivalent  water  system  since  the  internal  heat  transfer  coefficient 
is  now  known  for  a  wide  range  of  materials.  This  comparison,  made  for  a 
typical  advanced  interceptor  trajectory  (Figure  1),  shows  that  the  CONAP 
system  using  a  porous  tungsten  matrix  and  ammonia  coolant  is  50  percent 
lighter  than  a  water  system.  Figure  1  also  shows  that  the  CONAP  system, 
augmented  with  blind  holes,  is  40  percent  lighter  than  a  water  system. 

To  optimize  the  CONAP  system  efficiency  in  n  flight  environment,  the  heat 
transfer  coefficient  and  flow  characteristic  data  were  used  to  concep¬ 
tually  design  several  nose  tips  that  could  be  used  in  an  Advanced  Inter¬ 
ceptor  application.  The  mass  flow  variation  between  ideal  coolant 
requirements  and  actual  mass  flow  (Figure  2)  was  used  as  a  major  criterion 
for  selection  of  the  best  materials  and  concepts.  The  results  of  the  con¬ 
ceptual  design  activity  showed  that  use  of  the  Sylvania  and  Wah  Chang 
porous  tungsten  materials  (80  percent  by  volume)  resulted  in  the  lowest 
coolant  flow  variation.  Figure  '(  shows  a  flow  characteristic  correlation, 
the  reciprocal  of  the  permeability  versus  the  inertial  resistance  coeffi¬ 
cient,  of  all  of  the  materials  tested.  The  figure  identifies  those  materials 
which  fall  into  the  region  of  low  efficiency  due  to  the  excessive  sidewall 
mass  flow  that  could  not  be  effectively  reduced  even  with  a  maximum  of  20 
microns  of  tungsten  vapor  deposit.  The  region  of  high  material  efficiency 
also  has  limitations  due  to  very  high  chamber  pressure  requirements  (P  > 

6000  psi)  and  low  wall  thicknesses  (t  <  0.1  inches).  The  Sylvania  andC 
the  Wah  Chang  (80  percent  by  volume)  porous  tungsten  fall  within  the  opti¬ 
mum  range  shown  between  the  two  vertical  lines  in  Figure  3.  The  specific 
nose  tip  design  that  minimizes  the  mass  flow  variation  is  shown  in  Figure  4. 
This  design  consists  of  a  Wah  Chang  80  percent  porous  tungsten  material 
augmented  with  blind  holes  at  the  sonic  point  to  increase  the  coolant  flow, 
locally,  and  modified  with  a  20  micron  vapor  deposit  on  the  sidewall  to  de¬ 
crease  the  flow  rate,  locally. 

Scanning  electron  micrographs  (50X  to  10,000X)  of  the  materials  and  con¬ 
cepts  tested,  as  well  as  Schlieren  photographs  of  a  surrogate  flow  field 
exiting  the  tungsten  and  modified  tungsten  matricies,  basically  reinforce 
and  corroborate  the  data  correlations  and  conclusions  of  this  report. 

Based  on  the  results  of  this  evaluation,  a  thermal  test  program  consis¬ 
ting  of  plasma  arc  tests  over  a  range  of  heating  rates  and  stagnation  pres¬ 
sures  is  recommended  for  future  evaluation. 
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1 .  INTRODUCTION 


1 . 1  Background 

As  the  peak  velocity  of  an  ABM  interceptor  is  increased,  active  cooling 
systems  become  increasingly  attractive  for  use  in  nose  tips  and  air  vane 
leading  edges  (Reference  2).  At  currently  anticipated  burnout  velocities 
for  advanced  interceptors  the  combination  of  severe  thermal  and  pressure 
environments  produces  large  recession  rates  with  current  interceptor  nose 
tip  and  air  vane  leading  edge  materials.  Coupled  with  this  environment  is 
a  particle  environment  arising  from  both  natural  and  induced  causes  which 
further  increase  surface  erosion.  Large  recession  rates  directly  affect 
vehicle  performance.  Excessive  nose  tip  recession  will  create  a  vehicle 
center  of  pressure  shift  which  can  result  in  increased  control  system  weights. 
Air  vane  blunting  and  area  loss  cause  a  large  change  in  control  force  and 
hinge  moment  that  also  reflects  in  increased  control  system  weights. 

Candidate  systems  for  the  nose  tip  include  reinforced  phenolic  passive 
tips  and  transpiration  cooling.  Candidate  systems  for  the  air  vane  leading 
edge  include  passive  systems,  refractory  metals,  and  transpiiation  cooling, 
as  shown  in  Figure  1-1,  Transpiration  cooled  mettl  surfaces  offer  the 
advantages  of  reduced  erosion  because  of  their  inherently  high  resistance 
to  impact  damage. 

The  use  of  a  coolant  transpired  through  a  metal  matrix  wall  as  a  heat 
protection  surface  for  the  nose  tip  of  a  hypersonic  vehicle  is  not  new  and 
has  been  receiving  considerable  attention  in  recent  years  as  a  means  of 
achieving  a  shape-stable  nose  tip  for  high  ballistic  factor  reentry  vehicles 
(References  3,  A,  and  5).  Ground  and  flight  tests  have  been  conducted  on 
systems  with  water  as  the  coolant  and  stainless  steel  as  the  matrix  with 
the  coolant  introduced  to  the  wall  surface  through  either  discrete  holes  or 
by  distributed  porosity.  The  basic  principles  are  well  proven  and,  although 
there  are  uncertainties  in  the  downstream  cooling  effectivity  of  the  fluid 
film  layer  and  a  required  safety  factor,  transpiration  cooling  has  been 
used  in  practice  (Reference  5)  to  provide  zero-recession  nose  tips  under  the 
most  severe  ICBM  reentry  conditions. 

Although  water  has  a  relatively  high  efficiency  as  a  coolant  and  has  a 
good  storage  density,  the  work  of  Reference  6  shows  that  if  water  is  re¬ 
placed  as  a  coolant  with  a  reactive  gas  that  is  u^ed  with  a  high  surface 
temperature  matrix,  numerous  system  benefits  could  theoretically  be  realized. 
These  include  lighter  coolant  weights  and  lighter  expulsion  system  weights. 

A  material  that  can  operate  at  high  surface  temperature  is  an  ideal  matrix 
material  since  high  surface  temperatures  greatly  reduce  the  net  incoming 
heat  flux;  also,  with  a  material  that  can  operate  at  a  high  temperature  local 
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hot  spots  are  not  catastrophic.  A  porous  refractory  metal  such  as  tungsten 
would  be  an  excellent  material  for  this  application  except  that  at  high  sur¬ 
face  temperatures  tungsten  must  be  protected  from  attack  by  the  boundary 
layer  oxygen. 

The  CONAP  concept  which  combines  the  reactive  coolant,  ammonia,  with  the 
high  temperature  tungsten  matrix  depends  upon  the  ability  of  the  injected 
coolant  to  greatly  reduce  or  shut  off  the  flow  of  oxygen  reaching  the  re¬ 
fractory  hot  wall.  A  series  of  tests  were  conducted  to  demonstrate  experi¬ 
mentally  the  feasibility  of  this  concept  (Reference  7) .  The  results  of  these 
tests  showed  that  at  a  heat  flux  of  3500  Btu/ft2-s  and  an  enthalpy  of  4000 
Btu/lb,  nr-  measurable  recession  occurred  in  a  5-second  test  where  ammonia 
was  used  to  deplete  oxygen  at  the  surface  of  a  porous  tungsten  matrix.  As 
a  result  or  these  tests,  system  development  tasks  were  defined  to  provide 
the  data  required  to  perform  a  design  of  an  actual  tip  or  air  vane  leading 
edge.  These  tasks  included  the  evaluation  of  the  flow  characteristics  which 
control  the  mass  flow  distribution  within  the  porous  tungsten  matrix  and  the 
evaluation  of  the  internal  heat  transfer  coefficient  which  controls  the 
amount  of  heat  absorbed  by  the  coolant  as  it  passes  through  the  matrix.  The 
importance  of  knowing  the  internal  heat  transfer  coefficient  is  demonstrated 
in  Figure  1-2  which  summarizes  a  system  weight  comparison  between  a  water  sys¬ 
tem  and  a  CONAP  system.  The  comparison  is  made  parametrically  as  a  function 
of  gas  exit  temperature  which  has  a  direct  correspondence  to  the  internal  heat 
transfer  coefficient.  The  weight  comparison,  which  is  normalized  to  the  water 
system,  shows  that  if  heat  transfer  coefficients  are  high,  on  the  order  of 
1000  Btu/ft3-8-°R  at  the  flow  rates  of  interest  (around  5  lb/ft2-s),  a  weight 
savings  of  50  percent  could  be  realized  by  using  the  CONAP  system.  Low  inter¬ 
nal  heat  transfer  coefficients  of  CONAP,  however,  would  make  the  CONAP  concept 
less  attractive.  A  knowledge  of  the  values  of  heat  transfer  coefficient  as  a 
function  of  coolant  mass  flux  was  therefore  essential. 

Equal  importance  is  placed  on  the  definition  of  the  flow  characteristics  - 
permeability  and  inertial  resistance  coefficient  -  of  the  available  porous 
tungsten  materials.  The  flow  characteristics  of  a  material  and  the  modifica¬ 
tions  to  those  flow  characteristics,  such  as  blind  hole  augmentation  and 
surface  porosity  variations  through  vapor  deposition,  greatly  influence  the 
coolant  mass  flow  distribution  that  will  be  delivered. 

Once  these  parameters  are  defined  a  selection  of  the  optimum  materials 
can  be  made. 

1.2  Study  Objectives 

This  study  was  initiated  to  acquire  the  data  necessary  to  assess  the 
performance  of  a  nose  tip  or  air  vane  leading  edge  in  a  flight  environment. 

An  overall  experimental  approach  has  been  taken  with  a  series  of  laboratory 
tests  to  provide  the  required  data. 

The  specific  objectives  of  this  study  are  to: 

1  Determine  the  flow  characteristics  (permeability  and  inertial  resis¬ 
tance  coefficient)  of  all  materials  and  concepts  considered 
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2  Determine  the  internal  heat  transfer  coefficient  for  an  ammonia 
gas  passing  through  the  tungsten  matricies  considered 

3  Evaluate  and  optimize  several  tungsten  matrix  material  concepts 
which  include  a  wide  range  of  porous  materials,  materials  with 
discrete  holes,  materials  that  combine  porosity  with  discrete 
holes,  and  materials  whose  surface  is  vapor  deposited  with 
another  material 

4  Establish  CONAP  system  efficiency  for  these  material  concepts 
through  the  use  of  the  flow  characteristic  and  heat  transfer 
coefficient  data 

5  Select  those  concepts  and  those  combinations  of  concepts  which 
show  the  most  promise  for  future  system  use 

6  Define  the  test  model  configurations  for  further  work  in  plasma 

-  arc  tests  at  low  and  high  pressure. 

1 . 3  Approach 

The  test  plan  devised  to  meet  the  objectives  of  this  program  considers 
both  porous  tungsten  and  discrete  hole  tungsten  matricies.  Obtaining  and 
verifying  internal  heat  transfer  coefficients  and  flow  characteristics 
for  both  discrete  hole  and  porous  matricies  allows  the  greatest  latitude  of 
design  in  future  system  applications.  Discrete  hole  matrix  concepts, 
although  not  nearly  as  efficient  in  conserving  system  coolant  weights  as 
porous  materials,  have  much  greater  strength  capability.  In  the  event  that 
this  becomes  the  predominate  requirement  for  some  systems  due  to  its  anti¬ 
cipated  superior  erosion  capability,  data  were  taken  to  verify  the  fact 
that  standard  pipe  flow  and  heat  transfer  coefficients  can  be  used  or 
easily  modified  to  adequately  predict  both  ammonia  flow  coefficients  and 
heat  transfer  data  while  ammonia  is  flowing  through  a  discrete  hole  in  a 
hot  tungsten  matrix. 

In  the  flight  application  the  pressure  and  heat  transfer  distributions 
around  a  porous  nose  tip  and  air  vane  leading  edge  create  a  need  for  the 
coolant/matrix  system  to  deliver  only  a  required  amount  of  flow  to  each 
area.  Porous  materials,  although  very  efficient  as  compared  to  discrete 
holes,  could  be  even  more  efficient  if  the  flow  characteristics  could  be 
varied  within  the  matrix  to  minimize  the  mass  flow  variation  as  illustrated 
in  Figure  1-3.  No  materials  are  currently  available  with  a  variable  perme¬ 
ability  and  significant  development  is  required  to  create  this  capability. 
The  approach  taken  to  minimize  the  flow  variation  was  to  use  existing  mate¬ 
rials  and  modify  their  basic  flow  properties.  Consequently,  two  concepts 
were  added  to  the  basic  porous  matrix  concept  to  provide  a  more  efficient 
system.  Augmented  holes  were  considered  for  use  in  areas  where  the  flow 
requirements  were  high;  vapor  deposition  of  tungsten  was  considered  for  use 
in  areas  where  the  flow  rate  requirements  are  low.  Flow  augmentation  was 
accomplished  with  the  addition  of  blind  holes  which  were  not  drilled  com¬ 
pletely  through  the  material  but  only  drilled  partially  Into  the  matrix 
material.  This  concept  was  considered  because  it  was  anticipated  that  flow 
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rates  could  be  Increased  in  a  local  area  without  significant  degradation  of 
the  heat  transfer  coefficient  in  that  area.  The  second  method  by  which 
local  flow  can  be  controlled  is  by  depositing  tungsten  vapor  onto  the  sur¬ 
face  pores  and  decreasing  the  surface  porosity.  With  this  concept,  a  vapor 
is  deposited  on  the  surface  of  the  porous  tungsten  and  adheres  to  and  ef¬ 
fectively  enlarges  the  grains.  This  concept  decreases  the  flow  in  a  local 
area  without  greatly  affecting  the  heat  transfer  coefficient.  This  type 
of  i.oncept  is  useful  in  the  sidewall  region  of  a  nose  tip  or  air  vane  lead¬ 
ing  edge,  where  a  ae ’.tease  in  coolant  flow  is  desireable,  because  the  high 
pressure  potential  across  the  porous  wall  forces  more  coolant  through  the 
wall  than  is  required  to  neutralize  the  low  heat  rate.  Both  of  these  con¬ 
cepts  can  greatly  aid  in  tailoring  the  flow  to  minimize  the  variation  be¬ 
tween  actual  coolant  flow  ana  required  coolant  flow. 

Table  1-1  shows  the  material  concepts  on  which  heat  transfer  coefficient 
tests  and  the  flow  characteristics  tests  were  conducted. 

TABLE  1-1 
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2 .  MATERIAL  DESCRIPTION 


Two  types  of  commercially  available  porous  tungsten  materials  were  eval- 
l3ted  in  this  study.  Astro  Met  Associates  supplied  tungsten  material  that 
was  manufactured  to  final  shape  in  the  porous  state.  This  material  required 
no  further  processing  prior  to  use  as  test  specimen.  The  other  type  of  mat- 
erial  used  was  tungsten-copper ,  which  consisted  of  an  interconnected  pore 
tungsten  matrix,  infiltrated  with  copper.  This  material  was  received  in  bar 
stock  from  which  the  specimen  was  machined.  Further  processing  of  this  mat- 
crf^l  was  required  to  extract  the  copper j  however,  machining  was  performed 
prior  to  the  extraction  of  the  copper.  The  tungsten  copper  material  was  ob¬ 
tained  from  three  suppliers,  Wah  Chang,  Sylvania,  and  Mallory  Metallurgical. 


Copper  removal  was  performed  in  a  vacuum  furnace  in  which  the  samples 
we*e  inductively  heated  to  2550°F  in  a  0.5  x  10”5  mm  Hg  vacuum.  A  total  of 
four  hours  was  required  to  melt  and  evaporate  the  copper.  Prior  to  the  eva¬ 
poration  process,  and  at  intervals  during  the  process,  the  samples  were 
etched  with  a  nitric  acid-water  solution  to  expedite  the  copper  removal  pro¬ 
cess.  The  initial  etch  was  24  hours  while  subsequent  etches,  generally  three 
in  number,  were  about  1/4  hour  each.  A  detailed  description  of  the  equip¬ 
ment  used  can  be  found  in  Appendix  B. 

A  material  description  and  metallographic  evaluation  of  the  materials 
supplied  by  the  various  vendors  follows : 

2 . 1  Wah  Chang 

The  Wah  Chang  tungsten- copper  materials  obtained  were  80  volume  percent 
tungsten  —  20  volume  percent  copper,  and  70  volume  percent  tungsten  -  30 
volume  percent  copper.  Scanning  electron  micrographs  taken  of  the  surface 
of  these  Wah  Chang  materials  are  shown  in  Figures  2-1  through  2-6.  Figures 
2-1  through  2-3  show  the  80  percent  tungsten  augmented  with  blind  holes.  The 
tungsten  appears  very  uniform  in  structure  and  porosity,  and  no  evidence 
exists  that  the  placement  of  the  hole  in  the  sample  distorted  the  surface 
or  in  anyway  harmed  the  sample. 

Figures  2-4  through  2-6  are  similar  photographs  for  the  70  volume  per¬ 
cent  tungsten  material.  Again  the  uniformity  is  very  good,  although  not  as 
good  as  the  80  percent  tungsten.  For  both  materials,  the  pore  size  appears 
to  be  3  to  8  microns.  Although  these  samples  were  tested  at  high  temperatures, 
particle  surfaces  have  edges  which  are  quite  angular  and  appear  rough.  Ap¬ 
parently,  the  material  uniformity  permitted  adequate  gas  flow  to  prevent 
melting  or  smoothing  of  the  material  surfaces  by  diffusion. 
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Figure  2-1  .  Scanning  Electron  Micrograph  of  Surface 
of  Porous  Tungsten  Sample  Wah  Chang,  80  Volume 
Percent,  Showing  a  Typical  Blind  Hole;  6CX 


Figure  2-2.  Typical  Blind  Hole,  80  Percent 
Wah  Chang,  25 OX 


2 . 2  Sylvania 


Sylvania  also  produces  an  alloy  of  80  volume  percent  tungsten  -  20  volume 
percent  copper.  The  Sylvania  material.  Figures  2-7  and  2-8,  appears  as  uni¬ 
form  or  perhaps  even  more  uniform  than  the  Wah  Chang  tungsten.  The  surfaces 
are  rough  and  angular,  but  these  samples  were  not  heated  to  elevated  temper¬ 
atures.  Note  that  in  every  case  where  “untested”  is  specified,  elevated 
temperature  testing  was  not  done  on  the  particular  sample  photographed.  The 
original  particle  size  of  the  Sylvania  material  appears  to  be  much  the  same 
as  the  Wah  Chang  tungsten.  The  Sylvania  pore  size  seems  to  average  3  to  5 
microns  in  diameter. 

2.3  Astro  Met  Associates,  Inc. 

Astro  Met  makes  several  series  of  porous  materials  with  tungsten.  These 
series  (100,  300,  and  500)  represent  different  initial  powders  and  different 
proprietary  processing  methods.  Properties,  including  permeability  and  in¬ 
ertial  resistance  coefficient,  are  variable  within  ar.y  one  series. 

Three  types  of  tungsten  were  ordered  from  Astro  Met  in  the  form  of  cylin¬ 
ders,  1/2  inch  diameter  by  1/4  inch  thick.  Six  samples  of  each  type  of  series 
were  to  be  furnished,  final  machined  to  size.  The  types  were  10/15-60,  307- 
40,  and  307-60.  Although  processing  is  proprietary  to  Astro  Met,  the  follow¬ 
ing  information  is  known.  Type  10/15-60  is  a  100  series  material  with  15 
micron  pores  forming  interstital  tetrahedra  among  particles  roughly  blocky 
and  spherical  shape.  The  density  is  60  percent  of  theoretical.  This  mater¬ 
ial  was  fabricated,  machined  to  size,  and  received  for  test.  Type  307-40  is 
a  300  series  material  with  pores  about  7  mils  in  size,  forming  blocky  voids 
with  channels  about  0.1  the  diameter  of  the  voids.  The  particles  tend  to  be 
tetrahedral  in  shape  and  the  material  is  40  percent  of  theoretical  density. 

The  307-60  material  is  also  a  300  series  material  and  is  very  similar  to 
the  307-40  material.  The  density,  however,  is  60  percent  of  theoretical, 
the  permeability  is  lower,  and  the  inertial  resistance  coefficient  is  higher. 

Astro  Met  samples  are  shown  in  Figures  2-9  through  2-28.  Figures  2-9, 
2-10,  and  2-11  show  the  307-40  material.  In  addition  to  the  finer  pores  of 
10  to  20  microns,  shown  in  Figures  2-10  and  2-11,  pores  as  large  as  200  mic¬ 
rons  may  be  seen  in  Figures  2-9  and  2-10.  These  large  pores,  which  appear 
to  be  distributed  fairly  uniformly,  tend  to  increase  the  permeability  of  the 
specimen.  Figure  2-11  shows  the  tungsten  grain  boundaries,  which  do  not  nec¬ 
essarily  correspond  to  the  original  particle  boundary.  The  tungsten  surfaces 
are  relatively  smooth  and  rounded  in  comparison  to  previous  photographs  and, 
as  will  be  shown  later,  in  comparison  to  some  of  the  other  Astro  Met  material. 
This  implies  that  the  high  temperature  testing  may  have  caused  some  surface 
fusion  or  diffusion. 

The  10/15-60  material  is  shown  in  Figures  2-12,  2-13  and  2-14.  This  por¬ 
ous  tungsten  is  much  more  uniform  than  the  307-40  material  and  all  pores  seem 
to  be  approximately  15  microns  in  size.  The  surfaces  of  this  sample  are  also 
rounded  and  smooth.  In  Figure  2-14,  some  grain  boundaries  may  be  seen  al¬ 
though  generally  each  particle  is  one  grain.  Some  additional  grains  apparent¬ 
ly  resulted  from  recrystallization  at  the  points  where  particles  touched. 

The  flow  characteristics  of  this  material  allowed  excessive  coolant  flow 
rates  for  a  given  pressure  drop  across  the  specimen,  and  thus  the  material 
was  not  considered  for  ultimate  usage. 


The  J07-60  material  was  supplied  in  two  batches.  The  307-60  batch  No.  1 
material  is  shewn  in  Figures  2-15,  2-16  and  2-17.  Here  the  surface  is  smooth 
and  rounded.  This  porous  tungsten  is  similar  to  307-40  except  for  density. 
The  307-60  material  (Figure  2-15)  has  a  few  scattered  large  holes  of  a  max¬ 
imum  of  200  microns  in  diameter,  similar  to  the  large  pores  of  307-40  mater¬ 
ial.  Figure  2-15  shows  the  machining  marks  on  the  surface.  In  Figure  2-17, 
a  few  scattered  small  holes,  approximately  10  microns,  may  be  seen.  Again, 
the  structure  does  not  seem  to  be  uniform  and  homogeneous.  Grain  boundaries 
may  be  observed  in  Figure  2-17;  the  grains  appear  to  be  smaller  on  the  aver¬ 
age  than  the  particles  of  Figure  2-9.  The  flow  characteristics  of  the  307-40 
and  the  307-60  batch  No.  1  materials  are  comparable.  Figures  2-18  through 
2-20  illustrate  307-60  material  which  has  not  been  exposed  to  high  tempera¬ 
ture  testing.  Again,  the  surface  shows  a  rough  and  angular  nature.  Figures 
2-21  and  2-22  show  the  307-60  batch  No.  2  material  in  the  untested  condition. 
The  coarse  and  fine  holes  are  apparent  just  as  in  previous  samples.  The 
batch  No.  2  307-60  material  has  a  much  lower  permeability  and  a  much  higher 
inertial  resistance  coefficient  than  either  the  307-60  batch  No.  1  or  the 
307-40  materials,  as  will  be  shown  in  Chapter  3. 

To  show  the  effect  of  a  vapor  deposited  coating,  gold  was  applied  to  a 
sample  of  the  batch  No.  2  307*60  (Figures  2-23  and  2-24).  This  coating  has 
little  effect  on  tha  large  pores,  but  as  seen  in  Figure  2-24  the  coating 
fairly  well  plugs  any  of  the  very  fine  holes.  The  nodular -growth  nature  of 
the  vapor  deposit  is  apparent.  The  gold  coating,  as  expected,  decreased  the 
flow  through  the  pore  as  matrix. 

2.4  Mallory  Metallurgical  Company 

Mallory  produces  several  tungsten-copper  alloys;  eight  alloys  of  varying 
density  and  fabrication  techniques  were  supplied,  from  which  samples  were 
machined  and  the  copper  extracted  by  the  method  given  above.  Table  2-1  out¬ 
lines  the  sample  configurations.  Regular  sintering  was  done  at  temperatures 
less  than  240C°F,  whiie  high  temperature  sintering  was  done  at  about  2400°F. 

TABLE  2-1 

Tungsten-Copper  Material  Produced  by 
Mallory  Metallurgical 


Sample 

Designation 

Volume 

Percentage  of 
Tungsten 

Process 

Designation 

A 

r  50 

3W3  regular  temperature 

B 

50 

3W3  high  temperature 

C 

58 

10W3  high  temperature 

D 

58 

10W3  regular  temperature 

E 

64.6 

30W3  high  temperature 

F 

64.6 

30W3  regular  temperature 

G 

80 

50W3  regular  temperature 

H 

75 

40W3  high  temperature 
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Figure  2-15.  Astro  Met  307-60,  50X 


Figtve  2-16.  Astro  Met  307-60,  250X 


Figure  2-20.  Untested  Astro  Met  307-60,  1000X 


Figure  2-21.  Untested  Astro  Met  307-60 
Batch  2,  5 OX 


Figure  2-2?..  Untested  Astro  Met  307-60 
Batch  2,  1000X 


Figure  2-23.  Untested  Astro  Met  307-60,  Batch  2 
Vapor-Deposited  Gold  on  Surface,  50X 


Bate 


All  of  the  Mallory  material  is  very  uniform  in  nature.  Mallory  A,  Figures 
2-25  and  2-26,  appears  to  have  been  made  from  very  fine  powder  and  is  quite 
porous,  as  would  be  anticipated  from  the  low  density  of  50  percent.  The 
larger  pores  average  about  5  microns  in  size.  Mallory  C  material,  Figures 
2-27  and  2-28,  has  pores  of  about  the  same  size,  5  microns,  but  the  particle 
size  appears  larger,  apparently  resulting  from  sintering  at  the  high  tem¬ 
peratures  or  starting  with  a  coarser  particle.  In  Figure  2-27  some  surface 
contamination  is  seen.  Epoxy  was  used  to  bond  the  samples  into  a  holder  to 
take  Schlieren  photographs  and  removing  the  samples  afterwards  without  con¬ 
taminating  the  surface  was  difficult. 

Several  of  the  Mallory  materials  were  coated  from  1  to  5  microns  in  thick 
ness  by  tungsten  vapor  deposition.  The  Mallory  D  material  which  was  vapor 
deposited  with  tungsten  is  shown  in  Figures  2-29  through  2-32.  Again  the 
base  material  is  very  porous,  is  uniform,  and  has  pores  as  large  as  10  mi¬ 
crons.  The  nodular  nature  of  the  vapor  deposited  tungsten  is  readily  seen 
at  the  higher  magnifications,  particularly  in  Figures  2-31  and  2-32.  For 
comparison  the  uncoated  side  of  sample  D  is  shown  in  Figure  2-33,  which  may 
be  compared  directly  with  Figure  2-31 .  The  tungsten  vapor  deposit  reduced 
the  size  of  the  pores  in  the  primary  structure  by  enlarging  the  tungsten 
grains  themselves. 

Figures  2-34  and  2-35  illustrate  the  Mallory  E  material,  while  Figures 
2-36  and  2-37  show  the  H  tungsten.  Again,  some  contamination  is  present  in 
these  photographs.  Both  materials  are  uniform  with  pores  as  large  as  5  mi¬ 
crons.  As  with  most  Mallory  materials,  the  surface  of  the  particles  is 
fairly  smooth  and  regular,  whether  the  samples  were  heated  in  testing  or 
not.  The  G  material  is  similar  to  the  Wah  Chang  and  Sylvania  materials. 
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Figure  2-26 


Mallory  A,  1000X 


Figure  2-27.  Mallory 


Figure  2-28.  Mallory  C 


Figure  2-29.  Untested  Mallory  D,  1.5  Micron 
Vapor-Deposited  Tungsten  on 
Surface,  SOX 


On 


Figure  2-30.  Untested  Mallory  D,  1.5  Micron 
Vapor -Deposited  Tungsten  on 
Surface,  1000X 


Figure  2-31.  Untested  Mallory  D,  1.5  Micron 
Vapor -Deposited  Tungsten  on 
Surface,  5000X 


Figure  2-32.  Untested  Mallory  D,  1.5  Micron 
Vapor -Deposited  Tungsten  on 
Surface  5000X 


Figure  2-33.  Untested  Mallory  D,  5000X 


Figure  2-34 


Untested  Mallory  E,  50X 


Figure  2-35.  Untested  Mallory  E,  1000X 
Some  Surface  Contamination 


Figure  2-36.  Untested  Mallory  H,  5 OX 


Figure  2-37.  Untested  Mallory  H,  1000X 
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TEST  RESULTS 


3.1  Flow  Characteristics  Tests 

The  local  mass  flow  transpiring  through  a  porous  matrix  is  defined  by 
the  temperature  and  pressure  gradient  across  the  matrix  as  well  as  the  flow 
characteristics  of  the  matrix  itself.  While  the  pressure  and  temperature 
gradients  are  functions  of  the  flight  environment,  internal  heat  transfer 
coefficient,  coolant  thermophysical  properties,  and  inlet  pressure,  the 
matrix  permeability  (T)  and  inertial  resistance  coefficient  (6)  are  prop¬ 
erties  of  the  porous  matrix  alone.  Both  T  and  $  must  be  determined  experi¬ 
mentally  for  a  specific  material  and,  once  known,  define  the  flow  character¬ 
istics  of  the  material.  Since  T  and  6  are  properties  of  the  matrix  the 
flow  characteristics  tests  were  performed  using  nitrogen  over  a  pressure 
drop  range  of  100  to  800  psid  at  room  temperature.  A  description  of  the 
test  apparatus  and  test  procedure  and  a  discussion  of  the  data  reduction 
procedure  is  given  in  Appendix  A.  The  flow  characteristic  test  matrix  is 
shown  in  Table  3-1. 


TABLE  3-1 

Flow  Characteristics  Test  Conditions 


Concept 

Pressure 

Drop 

psig 

Matrix 

Thickness 

inches 

Sample 

1 

800 

1/4 

Porous  tungsten 

600 

400 

200 

2 

800 

1/4 

Porous  tungsten  augmented  with 

600 

blind  holes 

400 

200 

3 

600 

1/4 

Tungsten  copper  with  one 

discrete  hole 

4 

600 

1/4 

Tungsten  copper  with  several 

discrete  holes 

5 

800 

1/4 

Porous  tungsten  with  flow 

600 

control  by  vapor 

400 

deposition 

200 
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The  test  results  are  presented  in  Tables  3-II  through  3-V  referenced 

reSuitsa£«°thi“te  ug"T’  T“ble  3'H  Slm““lzes  the  flow  charecteristic 
?  unaltered  porous  tungsten  materials.  Values  of  6  and  T 

fit  to  the  Li f  re\Ult  fr°m  3  least  s<?uares»  straight  line  curve 
fit  to  the  cold  flow  data  obtained  on  each  specimen.  These  data,  consist¬ 
ing  of  mass  flow  rate,  gas  temperature,  and  pressure  drop  across  the  porous 
lunesten  specimen  are  combined  to  form  the  parameters  (P2-PQ2) /2piZRTX  and 
m/y.  When  this  first  parameter  is  plotted  as  the  ordinate  and  the  second  as 
the  abscissa,  the  reciprocal  of  permeability  is  the  intercept  of  a  straight 

line  Tabled  SVhJ*  ^  L*  resistance  coefficient  is  the  slope  of  the’ 
e.  Table  3  II  indicates  that  the  high  volume  density  materials  manu- 

by  Sylvaala  and  Wah  Chang  have  low  numerical  values  of  permeability 

that  i8h Jh  nU::  nrtia^  resistan:e  This  trend  Indicates 

that  a  high  pressure  drop  is  required  to  force  a  given  amount  of  flow  through 

coefficient  a  hlgher  value  of  Permeability  and  lower  resistance 

flowLh™ t  indicates  that  a  small  preset  drop  is  required  to  force  the 

rlZJ n  t  P°rOUS  matriX*  ThG  Mall°ry  serial  tested  had  the  widest 

the  oriLnarr  tCtC  ’  SinC6  the  V°W  Percenta8e  of  tungsten  in 

Iht  °r^ginal  tun8stan  c°PPer  ranged  from  50  to  75  percent.  Two  groups  of 

TouTall  th  2  We"  te8ted>  the  60  Percent  by  volume  tungsten  (307-60) 

group  and  the  40  percent  by  volume  tungsten  (307-40)  group.  The  307-60 

material  consisted  of  two  batches,  each  batch  representing  a  different 

manufacturing  process.  The  307-60  batch  No.  1  material  and  the  307-40 

h^tcb^1  S“bstanJially  have  the  same  resistance  to  flow,  although  the  307-60 
batch  No.  specimens  exhibited  slightly  lower  average  permeabilities  and 

rr^coe£flfe"t8 1** did  che  3o7-°  4~s«n..  s? fi« 1^- 

teristics  of  double  specimens  that  were  selected  from  each  group  were  also 
measured  and  are  tabulated  in  Table  3-II;  as  expected  little  change  in  6  or 
T  resulted.  The  flow  characteristics  of  the  Astro  Met  307-60  batch  No.  2 
material  are  vast  different  from  the  batch  No.  1  material.  Although  this 
shows  Astro  Met',  capability  to  vary  flow  characteristics  within  the 
material  group  by  changing  the  manufacturing  process,  the  data  reproduci¬ 
bility  from  batch  to  batch  has  not  been  verified  in  this  program. 

in  thrK°Ugh  3"v  Sh°W  SOme  °f  the  data  from  which  the  Parameters 

in  Table  3-II  have  been  taken.  Figure  3-1,  for  example,  shows  the  differ¬ 
ences  in  permeability  and  resistance  coefficient  of  two  samples  taken  from 
different  ends  of  an  80  percent  Sylvania  tungsten-copper  stock.  This  differ- 

standards^o^flight'hardware^se^8  t0  dCVel°P  material 

rid*  =sf . w^  «•- 

C  ang  70  percent  material  was  used  as  the  base  material  from  which  the  effect 

rid  £i™  sr 

a  e  shown  in  Table  3  III  and  Figure  3-6  respectively.  The  effect  on  flow 
characteristics  of  the  depth  of  the  blind  holes  is  tabulated  in  Table  3-III 
snd  shown  in  Figure  3-7.  Wah  Chang  80  percent  material  used  „  the  base 

holla*™  flow  Whjch.  ^he  e*fect  of  d£lPth  was  obtained.  The  effects  of  blind 
holes  on  flow  rate  through  a  porous  medium  can  best  be  seen  however  on  a 

curve  where  mass  flow  is  shown  as  a  function  of  pressure  drop  across  the 
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TABLE  3*11 


Flow  Characteristics  of  Porous  Tungsten  Materials  (Concept  1) 


Material 

Manufacturer 

Specimen 

Identification 

Percent 
Tungsten 
By  Volume 

Permeability 

(ft*) 

Resistance 

Coefficient 

(1/ft) 

Reference 

Figure 

Number 

Wah  Chang 

80 

1.37  x  10‘13 

1.93  x  108 

Sylvania 

End  1 

80 

1.26  x  10"13 

0.871  x  108 

3-1 

Sylvania 

End  2 

80 

0.89  x  10'13 

3.14  x  lo8 

3-1 

Wah  Chang 

End  1 

70 

3.27  x  10‘13 

3.02  x  107 

Wah  Chang 

End  2 

70 

0.61  x  10"13 

9.4  x  107 

Mallory 

A 

50 

0.36  x  lO"11 

0.21  x  106 

3-2 

Mallory 

B 

50 

0.14  x  lO’11 

0.52  x  106 

3-2 

Mallory 

C 

58 

0.19  x  lO'11 

1.04  x  106 

3-2 

Mallory 

D 

58 

0.65  x  10*12 

1.49  x  106 

3-2 

Mallory 

E 

65 

0.24  x  10  u 

0.47  x  106 

3-2 

Mallory 

F 

65 

0.25  x  10*12 

0.13  x  106 

3-2 

Mallory 

H 

75 

0.69  x  lo’13 

0.84  x  108 

3-2 

Astro  Met 

307-60(1)  Batch  2 

60 

0.17  x  10-13 

6.28  x  109 

3-3 

Astro  Met 

307-60(3)  Batch  2 

60 

0.63  x  10‘13 

1.33  x  109 

3-3 

Astro  Met 

307-60(4)  Batch  2 

60 

0.48  x  10'13 

1.54  x  109 

3-3 

Astro  Met 

307-60(6)  Batch  2 

60 

0.30  x  10'13 

1.42  x  109 

3-3 

Astro  Met 

307-60(7)  Batch  2 

60 

0.32  x  10‘13 

0.69  x  109 

3-3 

Astro  Met 

307-60(8)  Batch  2 

60 

0.26  x  10"13 

3.23  x  109 

3-4 

Astro  Met 

307-60(1)  Batch  1 

60 

0.458  x  10‘10 

0.62  x  106 

3-4 

Astro  Met 

307-60(2)  Batch  1 

60 

-10 

0.79  x  10  lu 

0.109  x  107 

3-4 

Astro  Met 

307-60(3)  Batch  1 

60 

0.111  X  10'9 

0.368  x  106 

3-4 

Astro  Met 

307-60(4)  Batch  1 

60 

0.165  x  10‘10 

0.299  x  '06 

3-4 

Astro  Met 

307-60(5)  Batch  1 

60 

0.364  x  10‘10 

0.432  x  106 

3-4 

Astro  Met 

307-60(6)  Batch  1 

60 

0.713  x  10'11 

0.195  x  106 

3-4 

Astro  Met 

307-60(7)  Batch  1 

60 

0.955  x  lO*11 

0.428  x  106 

3-5 

Astro  Met 

307-40(1) 

40 

0.158  x  10‘10 

0.507  x  106 

3-5 

Astro  Met 

307-40(1) 

40 

0.144  x  10“10 

0.224  x  106 

3-5 

Astro  Met 

307-40(2) 

40 

-10 

0.133  x  10 

0.195  x  106 

3-5 

Astro  Met 

307-40(3) 

40 

0.112  x  10‘10 

0.299  x  106 

3-5 

Astro  Met 

307-40(4) 

40 

0.115  x  10‘10 

0.365  x  106 

3-5 

Astro  Met 

307-40(5) 

40 

0.789  x  10*11 

0.242  x  106 

3-5 

Astro  Met 

307-40(5), 

40 

0.232  x  10*10 

0.139  x  106 

307-40(1)* 

■J  J 

Astro  Met 

307-40(6), 

40 

0.99  x  10"11 

0.184  x  106 

•XmK 

307-40(2)* 

Astro  Met 

307-60(1), 

60 

0.288  x  10"10 

0.425  x  106 

3-4 

307-60(5) 

Batch  1* 

Astro  Met 

307-60(2), 

60 

0.386  x  10'10 

0.534  x  106 

3-4 

307-60(6) 

Batch  1* 

*Double  Specimen 

3-3 
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TABLE  3-III 


Flow  Characteristics  of  Porous  Tungsten  Materials  Augmented  with 
Blind  Holes  (Concept  2) 


Percent 

Resistance 

Reference 

Material 

Tungsten 

Permeability 

Coefficient 

Figure 

Manufacturer 

By  Volume 

Hole  Configuration 

(ft2) 

(1/ft) 

Number 

Wah  Chang 

70 

1  blind  hole  0.020 
dia  x  0.15  inch  deep 

-13 

2.24  x  10 

2.84  x  106 

3-6 

Wah  Chang 

70 

5  blind  holes  0.020 
dia  x  0.15  inch  deep 

-13 

1.92  x  10 

-13 

1.61  x  10 

2.97  x  106 

1.86  x  106 

3-6 

Wah  Chang 

70 

9  blind  holes  0.020 
dia  x  0.15  inch  deep 

3-6 

-13 

1.48  x  10 

2.45  x  106 

Wah  Chang 

70 

17  blind  holes  0.020 
dia  x  0.15  inch  deep 

3-6 

Wah  Chang 

80 

25  blind  holes  0.020 
dia  x  0.030  in''!i  deep 

-13 

1.22  x  10 

2.96  x  107 

3-6 

Wah  Chang 

80 

25  blind  holer,  0.020 
dia  x  0.060  inch  deep 

-13 

1.47  x  10 

1.31  x  107 

3-7 

Wah  Chang 

80 

25  blind  holes  0.020 
dia  x  0.090  inch  deep 

-13 

1.83  x  10 

0.95  x  107 

3-7 

TABLE  3 -IV 


Effect  of  Flow  Direction  on  the  Flow  Characteristics  of  Porous 
Tungsten  Materials  Augmented  with  blind  Holes  (Concept  2) 
(Using  Wah  Chang  80  Percent  Tungsten  with  25  Blind  Holes) 


Hole 

Reference 

Depth 

Resistance 

Figure 

(inch) 

Flow  Direction 

Permeability 

Coefficient 

Number 

0.060 

Blind  holes  downstream 

-13 

1.47  x  10 

1.31  x  107 

3-8 

0.060 

Blind  holes  upstream 

-13 

1 .6  x  10 

2.03  x  107 

3-8 

0.090 

Blind  holes  downstream 

1.83  x  10*13 

0.95  x  107 

3-8 

0.090 

Blind  holes  upstream 

2.82  x  10"13 

1.47  x  107 

3-8 

TABLE  3-V 


Flow  Characteristics  of  Tungsten  with  Discrete  Holes 
(Concepts  3  and  4) 


Material 

Manufacturer 

Hole  Configuration 

Flow 

Coefficient 

Reynolds 

No. 

Wah  Chang 

5  through  holes 

0.025  inch  dia 

0.65 

5  x  103 

Wah  Chang 

5  cnrough  holes 

0.025  inch  dia 

0.70 

1.5  x  10* 

Wah  Chang 

5  through  holes 

0.030  inch  dia 

0.53 

5  x  10J 

Wah  Chang 

5  through  holes 
0.030  inch  dia 

0.55 

1.5  x  10* 

Wah  Chang 

1  through  hole 

0.030  inch  dia 

0.53 

5  x  10J 

Wah  Chang 

1  through  hole 

0.030  inch  dia 

0.55 

1.5  x  10* 

IIIKHW  HIMimi 
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m/y  x  10 

Figure  3-2.  Flow  Characteristics  of  Mallory  Tungsten 
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*■  80  percent  Wah  Chang 


Figure  3-7.  Flow  Characteristics  of  80  Percent  Wah  Chang  Tungsten  Augmented  with  Blind  Holes 


matrix.  Figures  3-8  and  3-9  show  this  effect.  Figure  3-8  shows  the  increase 
in  mass  flow  as  a  function  of  hole  depth.  For  example,  at  a  AP  *  600  paid, 
0.09  inch  holes  can  quadruple  the  flow  rate  in  80  percent  tungsten.  This 
technique  can  be  very  useful  if  a  large  amount  of  flow  is  required  in  a 
localized  area  because  of  very  high  heating  such  as  the  sonic  point. 

Figure  3-6  shows  the  effect  of  the  number  of  blind  holes  on  the  flow 
characteristics  of  porous  tungsten  matricies.  Flow  properties  of  the 
matricies  with  blind  holes  are  compared  with  an  identical  matrix  without 
blind  holes.  Specimens  with  1,  5,  9,  and  17  holes,  which  were  0.02  inch 
diameter  x  0.15  inch  deep,  were  drilled  in  a  0.25  inch  thick  porous  tungsten 
matrix.  The  specific  values  of  permeability  (the  reciprocal  of  the  inter¬ 
cept  on  the  ordinate  axis)  and  of  inertial  resistance  coefficient  (the  slope 
of  the  curve)  are  tabulated  in  Tables  3- II  and  3-III  for  each  case. 

The  presence  of  blind  holes  has  decreased  the  permeability  by  a  maximum 
of  50  percent  for  the  17  blind  hole  case  and  has  decreased  the  inertial 
resistance  coefficient  more  than  an  order  of  magnitude  for  all  cases.  Since 
the  slopes  of  the  curves  are  almost  the  same,  it  can  be  expected  that  at 
the  high  mass  flow  rates  of  interest  in  this  program,  m  >  0.5  lb/ft^-s, 
the  mass  flow-pressure  drop  relationship  will  show  a  very  weak  dependence 
on  the  number  of  blind  holes  in  a  matrix.  The  mass  flow-pressure  drop 
relationship  shown  in  Figure  3-9  verifies  this  observation.  It  is  apparent 
that  the  flow  rate  increase  resulting  from  the  addition  of  one  blind  hole 
is  so  large  that  additional  blind  holes  do  not  appreciably  enhance  the  flow 
any  further.  Although  a  detailed  flow  field  analysis  would  be  required  to 
fully  explain  this  phenomenon,  it  was  determined  that  this  analysis  was  not 
warranted  at  this  time. 

The  difference  in  mass  flow  for  a  given  pressure  drop  between  downstream 
blind  holes  and  upstream  blind  holes  is  shown  in  Figure  3-10.  The  flow 
characteristics  for  these  cases  are  tabulated  in  Table  3- IV.  There  is  a 
slight  difference  in  the  amount  of  flow  passing  through  a  specimen  at  a  given 
pressure  drop;  blind  holes  placed  on  the  upstream  side  are  more  effective. 
Because  of  geometrical  limitations,  upstream  blind  holes  can  only  be  easily 
placed  on  the  stagnation  point  of  a  nose  tip  or  on  a  stagnation  line  of  an 
air  vane,  leading  edge. 

With  the  data  from  Tables  3- III  and  3-IV,  the  effect  of  hole  depth  on 
the  inertial  resistance  coefficient  is  shown  in  Figure  3-11.  From  Figures 
3-7  and  3-11,  it  can  be  concluded  that  blind  hole  depth  has  little  effect 
on  permeability  but  has  a  large  effect  on  inertial  resistance  coefficient. 

The  characteristics  of  compressible  fluids  passing  through  pipes  are 
well  known  and  documented.  The  characterization  of  fluids  passing  through 
hot  tungsten  copper,  with  discrete  hole  diamettrs  of  interest  for  this 
study,  are  shown  in  Table  3-V  and  Figure  3-12.  The  flow,  as  expected,  can 
be  represented  in  terms  of  the  flow  function,  w/TR 

where: 

w  -  mass  flow  rate 

T  -  temperature 

R  -  gas  constant 

PT  -  pressure 

A  * 


flow  area 


Figure  3-8.  Effect  of  Blind  Hole  Depth  on  Flow  Through  80  Percent  Tungsten 


aP  -v  psfd 

Figure  3-9.  Effect  of  Blind  Hole  Number  on  Flow  Through  70  Percent  Tungsten 


Figure  3-10.  Effect  of  Upstream  and  Downstream  Blind  Holes  Flow 


1/ft 


Reynolds  Number,  NR  x  10-* 

Figure  3-12.  Flow  Characteristics  of  Discrete  Holes 
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Flow  coefficients  have  been  derived  for  the  0.025  inch  diameter  discrete 
holes  and  the  0.030  inch  diameter  holes  as  a  function  of  Reynolds  Number 
and  are  reported  in  Table  3-V.  The  flow  function  of  an  ideal  gas  with  a 
y  of  1.4  is  also  shown  in  Figure  3-12. 

The  data  presented  in  Table  3-VI  show  the  results  of  flow  control  by 
vapor  deposition.  The  purpose  of  this  concept  is  to  reduce  mass  flow  in  a 
region  in  which  it  is  not  needed.  For  example,  in  a  sidewall  region  of  a 
nose  tip  or  air  vane  leading  edge,  the  pressure  drop  is  very  high  because 
of  the  low  sidewall  local  pressure,  thereby  creating  the  potential  for  a 
high  sidewall  mass  flux.  However,  the  heat  rate  is  much  less  than  that 
found  in  the  stagnation  or  sonic  regions  and  a  high  mass  flow  rate  is  not 
necessary.  A  decrease  in  permeability  and  increase  of  resistance  coeffi¬ 
cient  can  be  achieved  through  vapor  deposition  which  in  effect  partially 
closes  the  surface  pores.  The  sputtering  technique  for  vapor  deposition 
is  described  in  Chapter  2  and  Appendix  B.  Both  tungsten  vapor  and  gold 
vapor  were  used  in  these  experiment. 

Although  gold  would  not  be  an  acceptable  material  for  the  flow  control 
of  a  CONAP  nose  tip  or  air  vane,  it  was  employed  for  these  cold  flow  ex¬ 
periments  in  order  to  determine  the  variation  in  flow  with  sputtering 
thickness  for  two  different  metals.  A  comparison  of  Figure  3-14  (tungsten) 
and  Figure  3-15  (gold)  at  a  thickness  of  1  micron  indicates  little  differ¬ 
ence  between  the  two. 


TABLE  3-VI 

Flow  Characteristics  of  Porous  Tungsten  Materials  with  Flow  Control 
By  Vapor  Deposition  (Concept  5) 


Material 

Manufacturer 

Specimen 

Identification 

Percent 

Tungsten 

By  Volume 

Vapor  Deposition 
Material  and 
Thickness 

Permeability 

Resistance 

Coefficient 

Mallory 

A 

50 

1  micron  tungsten 

0.22  x  10*11 

0.34  x  106 

Mallory 

C 

58 

2  micron  tungsten 

0.24  x  10"11 

3.62  x  106 

Mallory 

E 

65 

3  micfon  tungsten 

0.42  x  in”12 

2.25  x  106 

Mallory 

F 

65 

5  micron  tungsten 

0.23  x  10’12 

3.45  x  106 

Astro  Met 

307-60(6)  Batch  1 

60 

0.5  microu  gold 

2.03  x  lo'11 

3.85  x  105 

Astro  Met 

307-60(6)  Batch  1 

60 

1.0  micron  gold 

1.64  x  10'11 

5.55  x  105 

Astro  Met 

307-60(6)  Batch  1 

60 

1.5  micron  gold 

-12 

3.04  x  10  z 

2.89  x  105 

Mallory  A,  C,  E  and  F  materials  were  used  for  tungsten  vapor  deposition 
and  Astro  Met  307-60  batch  No.  1  material  was  used  for  gold  vapor  deposition. 
Flow  characteristic  data  is  tabulated  in  Table  3-VI.  Figure  3-13  shows  the 
flow  characteristics  of  the  Mallory  materials  after  vapor  deposition.  A 
comparison  of  mass  flow  rates  through  porous  tungsten  materials  with  and 
without  vapor  deposition  is  shown  in  Figures  3-14  and  3-15.  Figure  3-14 
shows  the  effect  of  tungsten  vapor  surface  deposits  on  the  Mallory  materials. 
Figure  3-15  shows  the  effect  of  progressive  0.5  micron  deposits  on  one  Astro 
Met  370-60  batch  No.  1  specimen.  For  example  at  200  paid  pressure  drop,  a 
1.5  micron  vapor  deposition  can  reduce  the  flow  rate  by  roughly  40  percent. 
Figures  3-16  and  3-17  show  how  permeability  and  resistance  coefficient  are 
affected  by  the  thickness  of  gold  vapor  deposition.  Note  that  the  “thick¬ 
ness”  is  measured  by  the  thickness  of  vapor  deposited  on  a  flat,  smooth 
tare  specimen  placed  next  to  the  work  piece. 
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Figure  3-14.  Effect  of  Tungsten  Coating  Thickness  on  Flow  through  Porous  Tungsten 


Pen.teabillty  x  Hf1*  -v  ft* 

Figure  3-16.  Effect  of  Gold  Coating  on  Permeability  of  Astromet  Porous  Tungsten 


1/ft 


The  porous  tungsten  material  flow  characteristics,  for  all  the  porous 
materials  tested,  can  be  represented  by  a  plot  of  the  reciprocal  of  the 
permeability  as  a  function  of  inertial  resistance  coefficient.  Although 
there  is  some  scatter,  the  data  is  consistent.  Figures  3-18  and  3-19  shows 
the  effect  of  the  volume  percentage  of  tungsten  on  the  inertial  resistance 
coefficient.  These  two  curves  indicate  a  median  value  of  expected  flow 
characteristics  properties  of  porous  tungsten  material  once  the  percent 
tungsten  by  volume  is  known. 

The  data  and  the  correlations  presented  in  this  part  of  the  report  can 
be  used  to  select  those  materials  whose  flow  characteristics  best  meet  the 
requirements  of  the  advanced  interceptor  mission. 

3.2.1  Heat  Transfer  Coefficient  Tests 

The  operating  temperature  of  a  CONAP  nose  tip  or  air  vane  leading  edge 
is  controlled  by  the  flight  environment,  the  ammonia  flow  rate,  and  the 
heat  transfer  rate  between  the  porous  tungsten  and  the  ammonia  coolant.  In 
the  ideal  case  of  thermodynamic  equilibrium,  the  heat  transfer  rate  between 
the  coolant  and  matrix  is  sufficiently  high  so  that  the  local  coolant  tem¬ 
perature  is  equal  to  the  local  matrix  temperature.  In  reality  this  does  not 
occur.  The  rate  of  heat  transfer  is  finite  and  is  controlled  by  the  internal 
heat  transfer  coefficient  (h1).  This  results  in  a  coolant  temperature  that 
is  lower  than  the  matrix  temperature  at  any  given  point  in  the  matrix.  The 
efficiency  of  the  CONAP  concept  depends  on  the  amount  of  heat  taken  from 
the  matrix  by  a  given  mass  of  coolant.  In  this  study,  laboratory  experi¬ 
ments  were  devised  to  obtain  the  data  necessary  to  calculate  the  internal 
heat  transfer  coefficient,  in  order  to  determine  the  effectiveness  of  the 
CONAP  concept.  A  test  fixture  was  built  whereby  an  ammonia  gas  stream  was 
passed  through  an  inductively  heated  porous  tungsten  plug.  Coolant  flow 
rates  ranged  from  0.5  to  8  lb/ft2-s  and  the  tungsten  plug  was  maintained 
at  temperatures  ranging  from  2000°R  to  5000*R  in  these  experiments.  The 
selection  of  these  test  conditions  was  based  on  estimates  of  required  cool¬ 
ant  flow  rates  and  matrix  temperature  conditions  on  the  nose  tip  of  an  ad¬ 
vanced  interceptor.  Measurements  taken  in  the  heat  transfer  experiment 
included  inlet  and  exit  coolant  temperature,  matrix  temperature,  coolant 
flow  rate,  and  upstream  pressure.  These  measured  data  were  then  factored 
into  a  computer  code  (Appendix  C)  which  computed  the  internal  heat  transfer 
coefficient.  A  detailed  discussion  of  the  test  description,  test  procedure, 
test  fixture  and  instrumentation  is  given  in  Appendix  A.  Test  specimens 
were  as  described  in  Table  3-1. 

The  heat  transfer  coefficient  data  for  the  porous  matrix  concepts  and 
the  discrete  hole  concepts  are  presented  separately.  Because  of  the  physi¬ 
cal  differences  in  these  two  cases,  the  data  in  each  had  to  be  correlated 
accordingly.  Heat  transfer  from  a  tungsten  matrix,  which  has  discrete 
holes,  to  a  coolant  gas  passing  through,  can  be  described  by  Newton's 
cooling  law  for  a  solid  body  convecting  heat  to  a  gas.  Heat  transfer  from 
a  porous  tungsten  matrix  to  a  coolant  gas  passing  through,  is  a  volumetric 
phenonemon  similar  to  heat  source.  Thus  the  surface  heat  transfer  phenomena 
of  the  discrete  hole  matrix  is  correlated  in  a  different  manner  from  the 
volumetric  heat  transfer  phenomena  taking  place  in  a  porous  matrix. 
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Figure  3-18.  Sintered  Tungsten  Flow  Correlation 
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Figure  3-19.  Effect  of  Volume  Percentage  on  Resistance  Coefficient 


For  the  case  of  the  porous  materials,  the  heat  transfer  coefficient 
that  was  computed  from  the  experimental  data  was  correlated  into  a  dimen¬ 
sionless  format.  The  data  were  put  into  the  form  of  Equation  (1). 


where: 

h1  =  heat  transfer  coefficient  (Btu/f t^-°R-s) 

L  ■  matrix  thickness  (feet) 

m  =  coolant  flow  rate  (lb/ft2-s) 

Cpc  =  coolant  specific  heat  (Btu/lb-°R) 

k  -  matrix  thermal  conductivity  (Btu/f t-s-°R) 
m 

A,B  *  curve  fit  coefficients 

For  each  specimen,  the  data  were  plotted  as  a  function  of 

h'L  ^  CPc 

- versus  — : - 

„  o  k 

Cn  m  m 

Pc 

and  a  least  squares  curve  fit  was  made  to  obtain  the  coefficients  in  Equa¬ 
tion  (1).  The  data  are  shown  in  Figure  3-20  through  3-28  and  Table  3-VII. 
Figures  3-20,  3-21,  and  3-22  show  the  data  and  the  curve  fit  for  the  straight 
porous  materials.  The  data  for  the  porous  specimens  augmented  with  blind 
holes  are  shown  in  Figures  3-23  through  3-26.  And  finally,  the  data  for  the 
porous  specimen  with  vapor  deposited  tungsten  coatings  are  shown  in  Figures 
3-27  and  3-28.  The  curve  fit  coefficients,  A  and  B,  as  shown  in  Equation 
(1),  corresponding  to  each  separated  specimen  are  tabulated  in  Table  3-VII. 

Table  3-VIII  tabulates  values  of  h*L/CpC  m  as  a  function  of  Lm  Cpc/Kn> 
as  shown  in  Equation  (2) ,  for  each  of  the  cu/ves  shown  in  Figures  3-20 
through  3-28.  Table  3-VIII  shows  the  range  of  heat  transfer  coefficient 
results  among  the  specimen  tested.  From  the  lowest  to  the  highest  value  of 
heat  transfer  coefficient,  the  variation  is  no  more  than  a  factor  of  three 
and  for  most  cases  it  is  much  less.  The  specimen  which  had  the  lowest  heat 
transfer  coefficient  at  all  values  of  Lm  Cp/1^  was  the  porous  tungsten  Wah 
Chang  specimen  augmented  with  one  blind  hole.  The  Mallory  A  specimen  had 
the  highest  heat  transfer  coefficient.  Table  3-VIII  also  shows  that  as  the 
number  of  blind  holes  are  increased  in  a  porous  tungsten  specimen, o the  heat 
transfer  coefficient  increases.  In  fact  at  the  higher  values  of  Lm  Cpc/km, 

which  correspond  to  flow  rates  ranging  from  4  to  8  lb/ft^-s,  the  heat  trans¬ 
fer  coefficients  for  the  nine  blind  hole  and  25  blind  hole  cases  were  com¬ 
parable  and  in  some  cases  higher  than  some  of  the  straight  porous  tungsten 
samples. 
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TABLE  3-VII 


Heat  Transfer  Coefficient  Correlations 


Material 

Manufacturer 

Specimen 

Identification 

Percent 

Tungsten 

By  Volume 

Blind  Hole 
Conf iguration 

Tungsten 

Vapor 

Deposition 

Thickness 

Coefficients* 

Reference 

Figure 

Number 

A 

B 

Mai  lory 

A 

50 

0.282 

0.0172 

3-20 

Mallory 

C 

58 

0.173 

0.0119 

3-21 

Wah  Chang 

- 

70 

0.262 

0.0063 

3-22 

Wah  Chang 

- 

70 

1  0.15  deep  x 

0.091 

0.0057 

3-23 

0,020  inch  dia 

Wah  Chang 

- 

70 

9  0.15  deep  x 

0.051 

0.029 

3-24 

0.020  inch  dia 

Wah  Chang 

- 

80 

25  0.06  deep  x 

O.COo 

0.037 

3-25 

0.020  inch  dia 

Wah  Chang 

- 

80 

25  0.09  deep  x 

0.071 

0.038 

3-26 

0.020  inch  dia 

Mallory 

A 

50 

1  micron 

0.308 

0.0 

3-27 

Mallory 

F 

65 

5  microns 

0.096 

0.014 

3-28 

*These  coefficients  fit  into  the  following  equation. 


TABLE  3-VIII 


Comparison  of  the  Heat  Transfer  Coefficients  of  the  Materials  Tested 

(h'L/CPc  m) 


Lm  CPc 

km 

Mallory 

A  Porous 

Mallory 

C  Porous 

Wah  Chang 
70  Porous 

Wah  Chang 

70  1  Blind 
Hole 

Wah  Chang 

70  9  Blind 
Holes 

Wah  Chang 

80  25  0.06 
Inch  Blind 
Holes 

Wah  Chang 

80  25  0.09 
Inch  Blind 
Holes 

Mallory 
A,  1u 
Coating 

Mallory 
F,  5u 
Coating 

2 

0.65 

0.39 

0.5 

0.22 

0.21 

0.31 

0.3 

0.6 

0.25 

4 

1.4 

0.89 

0.94 

0.46 

0.65 

0.9 

0.87 

1.25 

0.62 

6 

2.3 

1.45 

1.36 

0.73 

1.35 

1.8 

1.8 

1.85 

1.05 

8 

3.37 

2.15 

1.7 

1.1 

2.25 

3.0 

Extrapolated 

3.0 

2.46 

1.65 

TABLE  3 -IX 


Comparison  of  Test  Conditions  Under  Which  Heat  Transfer 
Coefficient  Data  Was  Taken 


Parameter 

Present  Study 

Grootenhuis 

Reference 

Schneider 
et  al 

Marco  and  Han 

Matrix  material 

Porous  tungsten 

Stainless 

steel 

Stainless 

steel 

Steel  wool 

Matrix  temperature 
range 

1000®R  -  5000 °R 

500 °R  - 
700°R 

700°R 

600  “R 

Coolant 

Gaseous  ammonia 
M  -  17 

Air 

M  -  29 

Water 

M  -  18 

Air 

M  -  29 

Mass  flow 
range 

0.5  -  8.0 
lbs/ft2-s 

7-16 
lbs/f t3-s 

1  -  8 
lbs/ft^-s 

0.1  -  0.6 
lbs/ft^-s 

Heat  transfer 

coefficient 

range 

2  -  2500 
Btu/ft3-»R-a 

600  -  1200 
Btu/ft3-°R-a 

550  -  1750 
3tu/ft3-s-*R 

0.2  -  3.5 
Btu/ft3-s-"R 
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k  comparison  of  the  data  presented  here  and  data  obtained  from  litera¬ 
ture  is  shown  in  Table  3-IX  and  Figure  3-29.  The  current  data  was  taken  at 
conditions  that  were  vastly  different  from  the  other  data.  The  data  taken 
from  literature  do  verify ,  however ,  that  the  trend  of  the  current  data  is 
correct.  Figure  3-29  shows  that  the  current  data  is  in  good  agreement  with 
that  of  Marco  and  Han,  Reference  8,  with  the  latter  being  at  the  low  end  of 
the  current  data.  The  current  data  however,  is  lower  than  Schneider,  Ref¬ 
erence  9,  who  used  water  as  a  coolant  in  a  stainless  steel  matrix.  This  is 
to  be  expected  since  liquids  have  higher  heat  transfer  coefficients  than 
gases.  Thus  the  Schneider  data,  the  Marco  and  Han  data,  and  the  current 
data  are  consistent.  The  data  of  Grootenhaus,  Reference  10,  who  used  hot 
air  that  was  cooled  by  flowing  through  a  stainless  steel  matrix,  agrees 
with  Schneider’s  data,  which  is  physically  not  consistent. 


The  discrete  hole  heat  transfer  data  are  compared  to  Nunner’s  equation 
(Reference  11),  for  heat  transfer  in  smooth  pipes.  This  equation. 


0.0396  (R  ) 

e 

d 


-1/4 


1  +  A(R  )'l/a  [Pr  (f/fj  -  1] 
ed 


(3) 


where : 


A  »  1.5  P 


-1/6 


shows  excellent  agreement  with  the  data  when  the  Prandtl  number  is  corrected 
for  ammonia  gas  as  shown  in  Figure  3-30.  The  data  itself  have  been  correlated 
as  a  function  of  Reynolds  number.  The  correlation  equation  used  is  shown  in 
Equation  (4) . 


Nu, 


A  R  +  B  R 


+  C  R 


(4) 


where : 


-  Nusselt  number  based  on  the  discrete  hole  diameter 

-  Reynolds  number  based  on  the  discrete  hole  diameter 

*  correlation  coefficients 

Table  3-X  gives  the  correlation  coefficients  that  fit  into  this  equation 
for  the  cases  tested. 

3.3  Schlieren  System  Test 

A  Schlieren  system  was  used  to  4 ‘view’ ’  the  flow  emerging  from  several 
of  the  configurational  concepts  that  have  been  tested  for  flow  character¬ 
istics  and  heat  transfer  coefficient.  The  Schlieren  system,  described  in 


Nu 


cd 

A,B,C 
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more  detail  in  Appendix  B,  is  a  system  whereby  density  gradients  in  the  flow 
field  are  distinguishable  and  can  be  viewed  on  a  screen  or  recorded  on  film. 
In  order  for  significant  density  gradients  at  room  temperature  to  exist  in 
subsonic  flow,  gases  with  a  wide  variation  of  molecular  weight  must  be  used. 
Thus,  helium  gas  was  used  as  the  coolant  in  the  tungsten  matrix.  The  helium 
gas  exited  the  matrix  into  air  at  room  temperature  ^...d  1  atm  pressure. 


TABLE  3-X 


Correlation  Coefficients  for  the  Heat  Transfer  to  Ammonia  Gas 
Passing  Through  A  Hot  Discrete  Hole  Tungsten  Matrix 


Discrete  Hole 

Correlation  Coefficients* 

Material 

Configuration 

A 

B 

C 

80  percent 
Wah  Chang 

5  through  holes 
0.025  inch  dia 

0.8137  x  10*2 

-2 

0.2945  x  10*6 

*6 

-0.1025  x  10'10 

-10 

80  percent 
Wah  Chang 

5  through  holes 
0.031  inch  dia 

0.99  x  10 

-2 

0.5239  x  10 

-6 

-0.2315  x  10 

-12 

80  percent 
Wah  Chang 

1  through  hole 
0.030  inch  dia 

1.031  x  10 

-0.4769  x  10 

0.3079  x  10 

*Correlation  coefficients  for  Equation  4. 


2 

Nu ,  ”  A  R  +BR  +CR 
d  Cd  Cd 


The  density  of  helium  is  roughly  one -seventh  of  air  at  the  same  pressure 
and  temperature,  thereby  giving  a  density  gradient  large  enough  for  the 
Schleiren  system  to  distinguish  flow  patterns.  The  computed  helium  mass 
flow  rates  are  much  lower  than  ammonia-nitrogen  flow  ra"es  expected  to  be 
used  in  a  flight  case  because  of  the  great  density  difference. 

All  of  the  materials  to  be  viewed  in  the  Schlieren  were  bonded  to  the 
holder  with  epoxy  around  the  periphery  of  the  specimen.  In  this  procedure 
emphasis  was  placed  on  the  prevention  of  leaks  between  the  holder  and  the 
edges  of  the  specimen.  Blank  specimen  were  used  to  assure  the  bonding 
process  had  in  fact,  sealed  the  sides  to  prevent  leaks  under  all  upstream 
pressure  conditions. 

Photos  were  taken  of  helium  through  an  open  pipe,  shown  in  Figure  3-31,  and 
photos  were  taken  of  helium  flow  through  a  discrete  hole  as  shown  in  Figures 
3-32  and  3-33.  These  runs  were  used  as  an  ‘  ‘index*  ’  for  the  other  cases  in 
this  experiment  in  the  sense  that  they  represent  limiting  cases  of  the  porous 
material  on  one  hand  and  the  blind  hole  concept  on  the  other.  Note  that  in 
Figure  3-33  which  shows  the  flow  through  a  matrix  with  5  through  holes  there 
is  obvious  coolant  starvation  in  the  areas  next  to  and  between  the  discrete 
holes.  Although  this  is  not  an  unexpected  result,  it  will  serve  as  a  contrast 
when  compared  to  the  porous  tungsten  matrix  augmented  with  blind  holes. 
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Figure  3-31. 


.  OPEN  PIPE 
2  *,■  f  =  5.6 

Speed  =  1/400  sec. 
m  =  0.3  x  10~4lb/sec 


/ 


OPEN  PIPE 
f  =  5.6 

Speed  =  1/400  sec. 
m  =  2.4  x  10~4  lb/sec 


Photograph,  Open  Pipe 


80  Percent 
Wah  Chang 
1  through  hole 
f  =  5.6 

Speed  =  1/400  sec. 
m  =  0.5  x  10-4  lb/sec. 


80  percent 
Wah  Chang 
1  through  hole 
f  =  5.6 

Speed  =  1/400  sec. 

A  =  2.4  x  10“3  lb/sec. 


Schlieren  Photograph,  80  Percent  Wah  Chang,  1  through  Hole 
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80  percent 
Wah  Chang 
5  thru  holes 
f  =  5.6 

Speed  =  1/400  sec. 
m  =  0.3  x  10~4  lb/sec. 


80-20  w-Cu 
Wah  Chang 
5  thru  noles 
f  =  5.6 

Speed  =  1/400  sec. 
m  =  1.2  x  10-4  ib/sec. 


Figure  3-33.  Schlieren  Photograph,  80  Percent  Wah  Ch.ing,  5  through  Holes 


The  flow  fields  emerging  from  the  Astro  Met  materials  are  shown  in 
Figures  3-34  through  3-37.  Astro  Met  307-40,  which  is  40  percent  by  volume 
tungsten  and  Astro  Met  307-60  batch  No.  1,  which  is  60  percent  by  volume 
tungsten,  not  only  have  the  same  flow  characteristics  but  also  seem  to 
have  downstream  flow  fields  that  are  similar.  It  appears  that  these  two 
materials  have  ‘‘discrete  holes’’  amidst  their  interconnected  porosity. 

This  is  borne  out  by  the  scanning  electron  microscope  photographs  which 
show  large  void  areas  which  probably  account  for  the  flow  behavior.  Flow 
through  the  Astro  Met  307-60  batch  No.  1  material  with  gold  coating  on  the 
downstream  side  seems  to  decrease  the  “discrete  hole’’  effect  of  the  earlier 
photos.  The  Astro  Met  307-60  batch  No.  2  material  appears  to  be  more  uniform 
than  the  batch  No.  2  material  in  the  sense  that  the  emerging  flow,  shown  in 
Figu-e  3-35,  does  not  have  “discrete  hole’’  effect  to  the  same  magnitude  as 
the  307-60  batch  No.  1  material.  This  is  not  unexpected  since  a  comparison 
of  Figures  3-23  and  3-25  shows  that  there  are  fewer  large  holes  in  the  Astro 
Met  307-60  batch  No.  2  material. 

Figure  4-38  shows  the  flow  emerging  from  the  Mallory  C  material.  Al¬ 
though  this  material  is  only  58  percent  tungsten  by  volume,  the  flow  is 
more  uniform  than  the  Astro  Met  materials.  Scanning  electron  microscope 
photographs  show  that  the  material  does  not  contain  a  significant  amount  of 
large  pores.  Figure  3-39  shows  the  flow  emerging  from  the  70  percent  Wah 
Chang  material.  The  uniform  flow  field  emerging  from  this  material  can  be 
expected  since  the  50X,  250X  and  1000X  photographs  in  Chapter  2  show  a  very 
uniform  material. 

Figures  3-40,  3-41,  and  3-42  show  the  flow  emerging  from  the  80  percent 
Wah  Chang  material  augmented  with  25  blind  holes,  the  70  percent  Wah  Chang 
material  augmented  with  9  blind  holes,  and  the  70  percent  Wah  Chang  material 
augmented  with  5  blind  holes.  In  all  cases  the  blind  holes  were  on  the 
downstream  (exit)  side  of  the  specimen.  These  figures  show  that  there  is  no 
flow  starvation  in  those  areas  where  there  are  no  holes.  This  is  an  obvious 
conclusion  when  these  photos  are  compared  with  the  discrete  flow  case  shown 
in  Figures  3-31,  3-32,  and  3-33  which  clearly  show  the  “no  flow  condition’’ 
between  the  discrete  holes.  The  fact  that  blind  holes  increase  the  local 
flows  without  starving  other  local  areas  can  be  used  to  great  advantage  in 
the  final  design  of  a  nose  tip  or  air  vane  leading  edge. 

A  porous  nose  tip  augmented  with  blind  holes  at  the  stagnation  point 
and  the  sonic  point  was  made  from  80  percent  Wah  Chang  material.  Helium 
flow  emerging  from  this  nose  tip  is  shown  in  Figure  3-43.  These  photographs 
show  that  there  is  flow  emerging  between  the  blind  holes  and  from  the  side- 
wall  . 
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Astro  Met  307-40(3) 
f  =  5.6 

Speed  =  1/400  s 
m  =  2.4  x  10-4  lb/sec 


Figure  3-34.  Schlieren  Photograph,  Astro  Met  307-40 
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Astro  Met  307-60(2) 
f  -  5.6 

Speed  =  1/400  s 
m  =  0.5  x  10~4  lb/sec 


Figure  3-35.  Schlieren  Photograph,  Astro  Met  307-60(2) 
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Figure  3-36. 


■Astro  Met  307-60(6) 

Batch  No .  1 
Gold  Coated 
f  =  5.6 

Speed  =  1/400  s 
Coating  Downstream 
m  =  0.3  x  10~4  lb/sec. 


Astro  Met  307-60(6) 
Gold  Coated 
f  =  5.6 

Speed  =  1/400  s 
Coating  Downstream 
m  =  2.4  x  10=4  lb/sec. 


Photograph,  Astro  Met  307-60,  Gold  Coated 


Figure  3-37.  Schlieren 


Astro  Met  307-60(1) 
Batch  No.  2 
f  =  5.6 

Speed  =  1/400  s 
m  =  2.4  x  10-4  lb/sec 


Astro  Met  307-60(1) 
Batch  No .  2 
f  =  5.6 

Speed  =  1/400  s 

rh  =  0.3  x  10“4  lb/sec 


Astro  Met  307-60(1) 


80  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s. 
m  =  0.3  x  10~4  lb/sec. 
25  Blind  Holes 


80  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s 
m  =  0.9  x  10”4  lb/sec. 
25  B.  ind  Holes 


Figure  3-40.  Schlieren  Photograph,  80  Percent  Wah  Chang,  25  Blind  Hole 


70  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s 
m  =  0.3  x  10-4  lb/sec. 
9  Blind  Holes 


70  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s 
m  =  1.1  x  10-4  lb/sec. 
9  Blind  Holes 


Figure  3-41.  Schlieren  Photograph,  70  Percent  Wah  Chang,  9  Blind  Holes 
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Figure  3-42.  Schlieren 


170  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s 
m  =  0.3  x  10-4  lb/sec. 
5  Blind  Holes 


170  percent  Wah  Chang 
f  =  5.6 

Speed  =  1/400  s 
m  =  1.2  x  10~4  lb/sec. 
5  Blind  Holes 


70  Percent  Wah  Chang,  5  Blind  Holes 


Figure  3-43. 


NOSE  CAP 

m  =  1.4  x  10~4  lb/se 


NOSE  CAP 

m  =  2.4  x  10-4  lb/se 


Photograph,  Nose  Cap 


4.  DESIGN  APPLICATIONS  AND  CONCLUSIONS 

The  basic  purpose  of  this  study  was  to  evaluate  transpiration  cooled 
materials  and  concepts  for  optimum  efficiency  in  a  flight  environment  as 
they  are  applied  to  a  nose  tip  and  an  air  vane  leading  edge.  The  results 
of  the  laboratory  testing  have  been  applied  to  a  design  of  a  nose  tip 
that  could  be  used  for  an  advanced  interceptor  flight  environment.  This 
analysis  has  shown  that  the  CONAP  concept  offers  a  weight  savings  over  a 
water  system.  In  addition  the  attempt  to1  minimize  the  mass  flow  variation 
between  ideal  coolant  requirements  and  actual  mass  flow  resulted  in  the 
definition  of  the  best  materials  and  design  concept.  The  Wah  Chang  and 
Sylvania  porous  tungsten  (80  percent  by  volume)  materials  were  selected. 
These  materials,  modified  by  blind  holes  at  the  sonic  point  and  vapor  de¬ 
position  on  the  sidewall,  were  selected  for  the  design  which  minimized  the 
flow  rate  variation. 


Several  conclusions  can  be  made  as  a  direct  result  of  the  data  taken. 

At  the  coolant  mass  flow  rates  expected  (5  lb/ft2-s)  for  a  typical  advanced 
interceptor  environment,  the  internal  heat  transfer  coefficient  is  on  the 
order  of  1000  Btu/f t3-s-°R.  The  impact  of  this  result  can  be  seen  in  Fig¬ 
ure  4-1  where  the  weight  of  a  CONAP  system  as  a  function  of  heat  transfer 
coefficient  is  compared  to  that  of  a  water  system.  For  this  case  the  CONAP 
concept  is  50  percent  lighter  than  the  water  system.  Figure  4-1  also  shows 
that  the  impact  of  the  blind  holes  on  system  weights  due  to  the  reduction 
of  internal  heat  transfer  coefficient  is  not  great. 

The  flow  characteristics  data  show  that  the  porous  matrices  augmented 
with  blind  holes  can  increase  the  coolant  mass  flow  rate  in  a  local  area. 

The  Schlieren  photographs  of  a  surrogate  coolant  showed  that  blind  holes  do 
not  starve  other  surface  areas  for  coolant.  The  combination  of  flow  char¬ 
acteristics  data  and  heat  transfer  coefficient  data  shows  that  the  process 
of  surface  vapor  deposition  decreases  coolant  mass  flow  in  a  local  area 
without  affecting  the  heat  transfer  coefficient.  These  modifications  to  a 
porous  tungsten  material  can  be  used  to  tailor  the  coolant  mass  flow  distri¬ 
bution  through  the  porous  matrix  to  better  meet  the  mass  flow  requirement 
distribution  of  a  nose  tip  or  air  vane  leading  edge.  The  blind  hole  technique 
can  be  used  in  the  sonic  point  region  of  a  nose  tip  or  stagnation  line  region 
of  an  air  vane  leading  edge  to  increase  the  local  mass  flow.  The  vapor  de¬ 
position  technique  can  be  used  on  the  sidewall  on  both  parts  to  decrease  the 
local  mass  flow. 

With  the  results  that  have  been  experimentally  obtained  it  is  now  possible 
to  define  the  best  of  the  available  tungsten  materials  that  can  be  used  in 
flight.  A  one  dimensional  analysis  was  conducted  using  the  results  of  the 
flow  characteristic  and  heat  transfer  coefficient  tests  to  determine  the 
most  nearly  optimum  material  and  concept.  A  typical  analysis  was  conducted 
in  the  following  manner.  First,  a  material  with  a  given  permeability  and 
inertial  resistance  coefficient  was  selected.  This  material  was  then  aug¬ 
mented  with  blind  holes  where  high  flow  rates  were  required,  and  tungsten 
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Normalized  System  Weight 


vapor  was  deposited  in  areas  where  the  flow  was  greater  than  desired.  A 
one -dimensional  coupled  heat  transfer  and  momentum  analysis  was  then  con¬ 
ducted  to  obtain  the  resultant  flow  rates  and  surface  temperatures.  These 
flow  rates  at  local  areas  then  were  plotted  as  a  function  of  the  non- 
dimensional  distance  from  the  stagnation  point.  In  addition  these  flow 
rates  can  be  compared  to  the  ideal  flow  rate  requirements  which  are  a  func¬ 
tion  of  the  heat  transfer  to  the  surface.  Figure  4-2  shows  the  flow  rate 
variation  between  the  ideal  coolant  flow  requirements  and  the  actual  coolant 
flow  that  will  be  delivered.  Table  4-1  summarizes  these  results.  Cases  1, 

2,  and  3  consider  modifications  (blind  holes  and  vapor  deposition)  to  those 
materials  that  can  be  readily  purchased  and  processed.  Of  these  the  most 
efficient  system  is  case  3,  since  the  flow  rate  variation  is  the  lowest. 

Figure  4-3  shows  this  conceptual  design  tnat  was  used  in  the  analysis.  This 
design  represents  the  minimum  variation  comparing  actual  with  desired  flow 
rates  that  could  be  theoretically  calculated  from  the  available  test  data 
for  available  materials. 

Figure  4-4  shows  the  flow  characteristic  comparison  of  all  of  the  porous 
tungsten  materials  tested.  The  figure  plots  the  reciprocal  of  the  perme¬ 
ability  versus  the  inertial  resistance  coefficient.  In  this  manner,  the 
materials  that  readily  allow  flow  through  them  appear  in  the  lower  left  hand 
corner  of  the  graph,  at  low  values  of  the  reciprocal  of  the  permeability  and 
low  values  of  the  inertial  resistance  coefficient  (B).  The  materials  which 
are  more  restrictive  to  the  flow  of  gases  through  them  appear  in  the  upper 
right  hand  corner  of  the  graph.  This  figure  shows  how  different  material 
flow  characteristics  can  affect  efficiency.  It  can  be  seen  that  low  ef¬ 
ficiency  can  result  with  a  material  that  allows  too  much  mass  flux  for  a 
given  pressure  drop  because  the  sidewall  flow  rates  will  be  much  greater 
than  required.  On  the  other  hand,  a  material  which  allows  too  little  flow 
for  a  given  pressure  drop  is  very  efficient  but  requires  very  high  chamber 
pressures  (>6000  psi)  and  very  thin  wall  thicknesses  (less  than  0.1  inch) 
to  deliver  the  flow.  Thus  the  optimum  flow  characteristic  region  for  mate¬ 
rials  is  in  between  these  other  regions  and  is  shown  in  Figure  4-4,  Two 
materials  tested  ha^e  flow  characteristics  which  fall  into  this  optimum 
region;  they  »iu  the  80  percent  Wah  Chang  and  the  80  percent  Sylvania  mat¬ 
erials  . 

The  results  cf  this  program  show  that  the  CONAP  concept  does  offer  po¬ 
tential  weight  savings  over  a  conventional  water  system.  The  next  phase  of 

this  development  would  logically  be  the  exposure  of  the  CONAP  concept  in  a 

plasma  arc  environment  over  a  range  of  pressures  and  enthalpies  that  are 
representative  of  currently  anticipated  advanced  interceptor  flight  regimes. 

In  the  current  program  a  one-dimensional  experimental  and  analytical  approach 
was  used  to  evaluate  and  select  the  most  nearly  optimum  materials  and  concepts. 
Testing  of  these  materials  and  concepts  in  the  plasma  arc  would  verify  cur¬ 
rent  conclusions  as  they  apply  to  three  dimensional  configurations.  It  is 
anticipated  that  the  three  dimensional  model  configuratoons  should  include 
a  discrete  hole  refractory  metal  alloy  model,  a  porous  tungsten  model,  a 
porous  tungsten  model  augmented  with  blind  holes  at  the  sonic  point  and  a 

porous  tungsten  model  augmented  with  blind  holes  at  the  sonic  point  and 

modified  with  vapor  deposition  of  tungsten  material  on  the  sidewall.  The 
plasma  arc  data  obtained  for  these  models  would  quantitatively  yield  the 
mass  flow  variation  between  ideal  requirements  and  actual  mass  flow.  This 
determines  the  mass  flow  safety  factor  requirements  for  each  case.  Cross- 
flow  effects,  the  CONAP  concept  performance  under  low  and  high  pressure 
environments  and  the  general  structural  integrity  of  the  model  are  also 
parameters  that  can  readily  be  examined  in  plasma  arc  testing. 
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RECIPROCAL  OF 
PERMEABILITY 

1/r  X  10UFT~2 


SYM  MATERIAL  DISC 

O  80%  WAH  CHANG 
□  80%  SYLVANIA  (1) 

O  >0  4  SYLVANIA  I2i 
^  75%  MALLORY  H 

A  70%  UUI'M  CHANG  (1) 

^  70%  WAH  CHANG  (2) 

O  58%  MALLORY  C 
A  58%  MALLORY  D 
n  50%  MALLORY  A 
O  50%  MALLORY  B 
X  (40%)  ASTROMET  307-40 
A  60%  ASTROMET  307-60 

BATCH  NO,  2  ^ 


MASS  FLOW 
MISMATCH  4 
TOO  GREAT 

LOW  EFFICIENCY 


BEST  OF  THE 
READILY  PURCHASED 
MATERIALS 


HIGH  EFFICIENCY 

CHAMBER 
PRESSURE 
REQUIREMENTS 
TOO  HIGH 
>600C  P'l 
AND  WALL 
THICKNESS  TOO 
LOW  <0.1  INCH 


100,000 


INERTIAL  RESISTANCE  COEFFICIENT  X  10*  FT1 

Figure  4-4.  Flow  Characteristic  Correlation  for  Several 
Porous  Tungsten  Materials 
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A. 1  Flow  Characteristics  Tests 


Rewriting  (5-4) ,  there  follows 


Since  all  values  shown  in  Equation  (5)  may  be  measured  at  each  test 
point,  with  the  exception  of  r  and  £,  we  may  plot  (p2  -  PQ2)/2  mZRTXy  (ordi¬ 
nate)  versus  m/y  (abscissa) .  Next  a  straight  line  may  be  fitted  to  the 
available  data  points.  The  straight  line  intercept  at  the  ordinate  axis 
is  equal  to  1 /T;  g  is  the  slope  of  the  straight  line. 

A  similar  method  is  employed  to  determine  r  and  3  if  a  liquid  test  cool¬ 
ant  is  used.  In  this  case,  however,  the  parameters  can  be  simplified  and  a 
p(P-PQ)/ymX  is  the  ordinate  ar.d  m/y  is  the  abscissa. 

A  Fortran  IV  computer  program  was  written  to  determine  I  and  g  employing 
these  relations.  The  Fortran  IV  listing  is  shown  in  Appendix  D. 

To  simultaneously  measure  the  quantities  required  to  obtain  r  and  g  for 
flow  through  several  specimens,  a  test  fixture  was  built.  A  schematic  dia¬ 
gram  of  the  test  fixture  is  shown  in  Figure  A-1 .  The  test  fixture  is 
designed  so  that  the  back  pressure  on  each  of  the  six  porous  tungsten  speci¬ 
mens  may  be  controlled  by  means  of  needle  valves  1  through  6.  Since  the 
upstream  pressure  on  all  six  specimens  is  equal,  each  pressure  drop  may  be 
set.  An  additional  circuit  is  also  available  on  each  of  the  six  specimens 
so  that  individual  flow  rates  and  discharge  pressures  may  be  measured.  For 
example,  if  the  coolant  flow  rate  and  specimen  discharge  pressure  through 
circuit  4  is  to  be  measured,  gate  valves  4A  and  4C  are  opened  and  4B  is 
closed.  This  will  force  the  coolant  to  pass  through  the  metering  venturi 
upstream  of  needle  valve  4.  During  fabrication  of  the  test  fixture,  care 
was  taken  to  equalize  the  two  circuits’  flow  pressure  drops.  The  flow 
metering  venturi  is  designed  so  that  most  of  the  pressure  drop  at  the 
throat  is  recovered  in  its  diffuser  section. 

The  test  fixture  fabricated  for  these  tests  is  shown  in  Figures  A-2  and 
A-3.  To  measure  the  required  test  parameters  it  was  necessary  to  include 
as  part  of  the  apparatus  a  mass  flow  measuring  device  and  pressure  measuring 
equipment.  The  following  paragraphs  briefly  describe  the  techniques  em¬ 
ployed  for  measuring  mass  flow  through  the  porous  tungsten  plugs. 

The  N2  supply  was  filtered  down  to  0.3  micron  to  prevent  clogging  of  the 
pores  in  the  tungsten  specimens.  Supply  pressure  levels  were  based  on  pre¬ 
dictions  of  expected  hot  wall  pressure  drops  across  0.25  and  0.50  inch  por¬ 
ous  walls  with  flow  characteristics  exhibited  by  the  specimens  received, 
and  a  downstream  pressure  of  6  atm.  The  resulting  range  of  pressures  varied 
from  400  to  800  psig  and  three  upstream  pressures  were  selected  for  testing: 
400,  600,  and  800  psig.  Since  flow  tests  were  to  be  performed  with  a  cold 
gas,  a  lower  pressure  drop  results  due  to  viscosity  effects. 
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The  schematic  shown  in  Figure  A-1  should  be  referenced  for  the  identity 
of  the  various  valves  discussed  in  the  detailed  test  procedure  below. 

The  cold  flow  test  procedure  is  as  follows: 

1  Install  a  porous  tungsten  specimen  of  the  type  shown  in  Figure  A-1 

~  of  each  of  the  six  holders. 

2  Set  supply  manifold  pressure  to  400  psig; 

3  Slowly  open  valves  1A  then  1C; 

4  Open  needle  valve  1  until  pressure  downstream  of  the  porous  plug 

-  is  360  psig; 

5  Open  valve  IB  and  close  1A  and  1C; 

6  Open  valve  6A  then  6C; 

7  Open  needle  valve  6  to  the  fully  opened  position; 

8  Based  on  the  downstream  pressure  reading  in  circuit  6,  compute 

-  desired  downstream  pressure  for  circuits  2,  3,  4,  and  5  by  deter¬ 
mining  a  uniform  pressure  spread  for  each  of  the  four  circuits; 

9  Open  valve  6B  and  close  6A  and  6C; 

10  Open  valve  2A  then  2C; 

11  Adjust  needle  valve  2  until  the  desired  downstream  pressure 
computed  in  step  7  is  achieved; 

12  Open  valve  2B  and  close  2A  and  2C; 

13  Repeat  steps  10  through  12  for  each  of  the  remaining  circuits  3, 

4 ,  and  5  . 

During  each  test,  some  readjustment  of  manifold  pressure  was  required  to 
maintain  the  desired  pressure  of  400  psig.  Once  the  downstream  pressure  was 
achieved  on  each  of  the  six  circuits,  the  data  was  recorded.  This  procedure 
was  repeated  for  the  600  psig  and  800  psig  manifold  pressure  levels.  The 
same  procedure  was  used  for  the  variable  thickness  and  permeability  tests. 

To  ensure  that  gas  leakage  by  the  specimens  did  not  occur  in  the  holders, 
a  solid  specimen  was  installed  in  each  holder  and  the  system  pressure  up¬ 
stream  of  the  disc  increased  to  approximately  1500  psi.  The  nitrogen  supply 
was  then  shut  off  from  the  test  apparatus  and  leakage  deduced  from  the  pres¬ 
sure  bleed  down.  Holder  modifications  were  made  until  there  was  no  decay 
in  the  pressure  level. 
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To  determine  the  characteristics  of  the  porous  tungsten  specimens  it  was 
necessary  to  determine  mass  flow  rate,  pressure  drop  across  the  specimens, 
flow  area,  nitrogen  temperature,  and  specimen  thickness.  These  data  were 
determined  as  described  below. 

Porous  tungsten  upstream  pressure  was  determined  employing  the  Heise 
gage  shown  in  Figure  A-2.  Downstream  pressure  was  measured  by  a  0  to  1000 
psi  pressure  transducer  in  conjunction  with  a  digital  millivolt  meter  shown 
in  Figure  A-3 .  The  transducer  output  was  also  recorded  on  a  millivolt  re¬ 
corder.  Transducer  calibration  curves  are  shown  in  Figure  A-4  and  A-5. 

Nitrogen  mass  flow  rate  was  determined  by  either  of  two  calibrated 
venturies.  For  purposes  of  these  tests  each  of  the  venturies  was  calib¬ 
rated  using  nitrogen  and  employing  flow  meters  of  known  accuracy.  Differ¬ 
ential  pressure  across  the  venturies  was  measured  with  either  of  two 
differential  pressure  transducers  depending  on  the  flow  rate  expected. 

Calibration  curves  for  these  transducers  are  shown  in  Figure  A-6  and 
A-7.  Mass  flow  was  computed  by  determining  theoretical  flow  as  a  function 
of  Mach  number,  calculating  ideal  Reynolds  number  based  on  the  throat  dia¬ 
meter  01  the  venturi,  finding  flow  coefficient  (Cp)  for  the  specific  venturi 
from  Figure  A-8,  and  finally  calculating  flow  as  the  product  of  Cp  and  the 
theoretical  mass  flow.  Note  that  the  value  of  CF  for  the  large  venturi  in 
Figure  A-8  is  greater  than  one.  This  behavior  is  caused  by  the  slightly 
elliptical  cross  section  of  the  throat  on  this  venturi.  Since  the  throat 
diameter  is  determined  by  measurement  with  a  pen  gage,  the  resultant  area 
computed  is  less  than  the  true  area.  This  results  in  a  low  theoretical  flow 
rate  thereby  making  the  Cp  correction  greater  than  one. 

A. 2  Heat  Transfer  Coefficient  "csts 

The  internal  heat  transfer  coefficient  tests  were  conducted  using 
ammonia  a^  the  coolant.  The  overall  test  schematic  diagram  is  shown  in 
Figure  A-9.  Coolant  is  supplied  from  the  coolant  supply  system  to  the 
test  fixture.  Upon  passing  through  the  hot  porous  plug  the  coolant  vents 
into  the  fume  hood  where  it  mixes  with  a  purge  supply  of  nitrogen.  From 
the  fume  hood  the  gaseous  mixture  is  passed  through  a  water  tank  where  the 
ammonia  is  absorbed  by  the  water  and  flushed. 

A  schematic  of  the  test  fixture  in  which  the  hot  porous  plug  is  heated 
and  through  which  the  coolant  passes  is  shown  in  Figure  A-10.  Coolant  was 
supplied  through  a  flexible  hose,  passed  through  a  quartz  phenolic  tube, 
and  into  the  boron  nitride  specimen  holder.  The  boron  nitride  specimen 
holder  is  capable  of  withstanding  the  5000°R  tungsten  temperature  and  at 
the  same  time  is  not  electrically  coupled  to  the  RF  heater. 

The  boron  nitride  specimen  holder  is  attached  to  the  quartz  phenolic 
inlet  tube  by  means  of  threaded  quartz  phenolic  retainer  rings.  The  system 
has  an  “0”  ring  and  grafoil  washers  to  prevent  leakage.  Figure  A-11  shows 
the  detailed  drawings  of  the  test  fixture  hardware.  Figure  A-12  shows  the 
actual  hardware  as  assembled  and  in  an  exploded  view.  Figure  A-13  shows  the 
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A* 4.  Downstream  Pressure  Transducer  Calibration  Curve 
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Figure  A-7.  Differential  Pressure  Transducer  Calibration  (High  Range) 
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Figure  A-9 .  Heat  Transfer  Coefficient  Test  Schematic  Diagram 
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laboratory  placement  of:  the  components  of  the  heat  transfer  coefficient 
tests.  These  components  include  the  induction  heater,  optical  pyrometer, 
coolant  supply  system,  and  the  millivolt  and  chart  recorders.  Figure  A-13B 
shows  the  placement  of  the  heat  transfer  coefficient  test  assembly  inside 
the  heater  coils.  The  assembly  was  housed  in  a  transite  box  which  is  in¬ 
side  the  fume  hood.  Also  seen  in  Figure  A-13B  is  the  thermocouple  mechanism 
used  to  insert  the  thermocouple  into  the  exit  of  the  boron  nitride  holder. 

Thermocouples  were  employed  to  measure  coolant  inlet  and  exit  tempera¬ 
tures.  An  optical  pyrometer,  as  shown  in  Figure  A- 10,  was  used  to  measure 
tungsten  temperature  at  the  exit.  Coolant  A-10  flow  rate  and  inlet  pres¬ 
sure  level  were  determined  employing  the  calibrated  venturi  and  the  pressure 
transducers  that  are  part  of  the  coolant  expulsion  system.  Three  flow  mea¬ 
suring  venturies  were  used  for  flow  measurement.  Tabel  A-I  indicates 
the  various  venturies  and  differential  pressure  levels  expected  during  the 
tests.  Once  the  test  data  discussed  above  has  been  obtained,  the  internal 
heat  transfer  rate  between  the  coolant  can  be  determined  employing  the  analyt¬ 
ical  technique  explained  in  Appendix  C.  The  porous  tungsten’s  permeability 
(F)  and  inertial  resistance  coefficient  (B)  determine  the  upstream  pressure 
necessary  to  force  the  desired  coolant  flow  through  the  porous  tungsten. 

This  analytical  technique  has  been  modified  slightly  so  that  cold  flow  pres¬ 
sure  drop  and  alternatively  F  and  6  may  be  determined  if  pressure  drop  is 
known. 

TABLE  A-I 


Flow  Venturi  Selection 


Coolant  Flow  Rate 
(lb/f t^-s) 

0.5 

1 

2 

5 

8 

Venturi  req  for  NH3 

Small 

Small 

Med 

Med 

Differential  pressure 
(psid) 

m 

0.5 

26 

18 

40 

Venturi  req  for  N2 

Small 

Med 

Med 

Large 

Large 

Differential  pressure 
(psid) 

13 

7 

22 

1 

2 

In  order  that  an  accurate  assessment  of  the  value  of  the  internal  heat 
transfer  coefficient,  h’ ,  could  be  made,  emphasis  was  placed  on  the  preven¬ 
tion  of  coolant  leakage  past  the  holder -specimen  interface.  To  assess  the 
leakage  past  the  grafoil  seals  placed  on  each  side  of  the  specimen,  both 
cold  flow  and  hot  flow  tests  were  run.  A  solid  tantalum  disk  was  fabricated 
and  used  in  the  tests.  The  test  procedure  consisted  of  assembling  the  hold¬ 
er  using  the  solid  tantalum  disk  and  pressurizing  the  system  and  holder  up¬ 
stream  of  the  disk  to  approximately  1000  psia.  A  gate  valve  just  upstream 
of  the  holder  was  closed  and  a  pressure  gage  attached  at  the  holder  was 
monitored.  No  decrease  in  pressure  was  observed  during  the  15  minute  test 
duration.  Next,  the  RF  generator  was  energized  and  the  tantalum  disk  heated 
to  approximately  5000°F  and  the  pressure  inside  the  holder  was  again  observ¬ 
ed  for  approximately  15  minutes.  No  pressure  drop  resulted. 
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In  order  to  measure  coolant  discharge  temperature,  a  traversing, 
radiation-shielded  thermocouple  was  fabricated  and  installed.  Thi.s 
thermocouple  was  platinum,  platinum  -  13  percent  rhodium,  capable  of  meas¬ 
uring  temperatures  to  3100°F.  The  influence  of  both  the  RF  field  used  for 
specimen  heating  and  the  porous  specimen  radiation  effects  on  the  thermo¬ 
couple’s  measurement  of  coolant  discharge  temperature  were  examined. 

Since  the  thermocouple  itself  is  metal,  it  can  couple  to  the  RF  field 
and  be  heated  inductively.  In  addition,  the  high  field  voltage  could  also 
influence  the  thermocouple  measurements.  These  effects  were  eliminated  by 
designing  the  specimen  holder  so  that  the  thermocouple  was  outside  the  field 
of  the  inductive  coil.  Two  tests  were  conducted  to  determine  if  thermo¬ 
couple  coupling  had  been  eliminated.  These  are  described  below: 

1  The  test  facility  was  assembled  without  a  porous  specimen  in  the 
~  holder.  The  traversing  thermocouple  was  located  in  its  position 

for  measurement  of  gas  discharge  temperature.  The  inductive  heater 
was  turned  on  and  the  work  coil  actuated  to  the  maximum  current 
(450  amps).  Since  no  coolant  was  being  heated,  all  increase  in 
thermocouple  temperature  would  be  due  to  RF  field  effects  on  the 
thermocouple  itself.  No  increase  in  temperature  was  measured. 

2  A  second  test  was  conducted  to  evaluate  accuracy  of  the  traversing 
“  thermocouple.  During  several  hot  flow  tests  the  coolant  discharge 

thermocouple  was  left  in  position  and  the  RF  field  switched  off. 

No  sudden  change  in  coolant  discharge  temperature  occurred.  If 
there  had  been  coupling  between  the  RF  field  and  thermocouple,  a 
sudden  decrease  in  indicated  coolant  temperature  would  have  resulted. 

The  influence  of  specimen  radiation  on  traversing  thermocouple  measure¬ 
ment  was  also  evaluated.  The  following  tests  were  run: 

1  A  special  holder  was  designed  with  a  90  degree  bend  in  the  coolant 

~  exit  duct.  A  second  thermocouple  was  installed  in  the  holder  down¬ 

stream  of  the  90  degree  bend  so  that  it  could  not  be  affected  by 
specimen  radiation.  A  hot  flow  test  was  run  and  the  two  coolant 
temperature  measurements  compared.  No  difference  between  the  two 
was  apparent. 

2  A  second  test  consisted  of  traversing  the  thermocouple  into  position 
-  behind  a  solid  specimen  that  was  being  heated  inductively.  Since 

there  was  no  coolant  flow,  all  temperature  rise  sensed  by  the 
thermocouple  would  be  by  radiation.  No  temperature  increase  was 
measured . 

Tungsten  surface  temperature  was  measured  by  means  of  a  Leeds  and 
Northrup  automatic  optical  pyrometer,  which  operates  by  a  continuous  auto¬ 
matic  matching  of  the  radiation  intensity  of  the  target  with  a  standard 
variable  output  lamp.  Measurement  of  the  voltage  to  this  standard  lamp 
provides  a  direct  indication  of  the  brightness  of  the  target  and  the  resulting 
brightness  temperature.  The  calibration  tables  employed  for  converting  the 
millivolt  input  to  the  standard  lamp  to  brightness  temperature  were  supplied 
by  Leeds  and  Northrup. 
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In  orde*  to  convert  brightness  temperature  to  true  temperature  it  was 
necessary  to  consider  the  normal  spectral  emittance  of  the  target  material 
ac  0.65  micron.  These  data  were  obtained  from  Reference  12  and  employed  in 
the  following  relation  to  convert  brightness  temperature  to  thermodynamic 
temperature. 


where: 

Tg  =  thermodynamic  temperature  of  the  surface 

TB  =  bri8htness  temperature  of  the  surface 

A  =  wavelength  of  the  optical  pyrometer 

C2  =  second  radiation  constant  (2.5896  x  10^  °Ry) 

=  emittance  of  surface  at  A. 

The  automatic  optical  pyrometer  views  the  porous  tungsten  test  specimens 
through  a  pyrex  window.  Brightness  temperature  was  corrected  for  the  effect 
of  this  window  based  on  Figure  108  from  Reference  13. 

During  operation  the  test  chamber  was  filled  with  N2  or  a  mixture  of  No 
and  NH3.  References  14,  15,  and  16  indicated  that  there  is  no  effect  of  the 
N2  or  NH3  vapors  on  the  optical  temperature  measurement  of  tungsten. 
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B , 1  30  kW  Induction  Heater 

Westinghouse  Industrial  Radio  Frequency  Generators  are  specifically 
designed  for  industrial  induction  heating  applications.  They  conform  to 
NEMA  Standards  for  induction  and  dielectric  equipment,  AIEE  Proposed 
Standard  and  Test  Code  for  Induction  and  Dielectric  Heating,  applicable 
portions  of  the  National  Electric  Code,  and  with  current  JIC  Standards  as 
interpreted  by  I.H.E.A. ,  and  to  meet  the  maximum  radiation  requirements 
of  Part  18  of  FCC  Rules  and  Regulations.  The  generator,  designed  to  operate 
from  a  nominal  240  or  480  volt,  3”phase  system,  has  the  following  character 
istics: 


RATED  OUTPUT 

Power  (into  NEMA  Load) 

Frequency 
Amperes 
Voltage  (Max) 

INPUT  (AT  RATED  OUTPUT) 

kVA 

Frequency 

Phase 

Veits 

Po*»°r  factor 

OPERATING  CONDITIONS 

Ambient  temperature 
Cooling  water 
pressure  (min) 

Flow 

Inlet  temperature  (max) 

Inlet  pressure  (max) 

Purity  "  (Spec  T.ES  1000) 
Recommended  auxiliary  fuses 
For  440  volt  operation 
For  220  volt  operation 
Recommended  control  circuit  fuse 
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30  kW  30  kW 

450  kHz  4.0  MHz 

168  Amps  350  Amps 
4900  Volts  9000  Volts 


66  kVA 
60  Cycles 
3  Phase 
220/440 
0.95 


104°F  (40°C) 

30  psi  diff 

17  Gpm  at  30  psi  diff 

90°F 

90  psi 

IES  1000 

10A 

20A 

15A 


— . . 
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B.2  Coolant  Supply  System 


A  coolant  supply  system  compatible  with  both  nitrogen  and  anhydrous 
ammonia  is  shown  schematically  in  Figure  B-1.  Figure  B-2  shows  the  supply 
system.  This  system  consists  of  a  high-pressure  accumulator  employing  a 
free  piston  to  pressurize  the  liquid  ammonia,  a  solenoid  shut-off  valve,  a 
two-stage  filtering  system,  and  adjustable  needle  valve,  and  a  flow-measuring 
venturi.  During  the  nitrogen  tests,  the  accumulator  can  be  bypassed  and  the 
high-pressure  nitrogen  supplied  directly  upstream  of  the  solenoid  shut-off 
valve.  Prior  to  a  test,  the  metering  valve  downstream  of  the  final  filter 
is  preset  to  produce  the  desired  coolant  flow  under  room- temperature  condi¬ 
tions.  Expected  changes  in  coolant  viscosity  and  density  during  the  hot 
tests  are  accounted  for  in  calculation  of  desired  coolant  flow  setting  at 
room  temperature.  The  flow-metering  valve  is  located  as  close  as  possible 
to  the  model  holder  to  minimize  coolant  flow  response  time. 

B.3  Vacuum  Furnace 

A  vacuum  furnace  was  employed  to  boil  off  the  copper  from  the  porous 
tungsten  specimens.  This  furnace  (shown  in  Figure  B-3)  consists  of  a  water 
cooled  stainless  steel  tank  with  a  viewing  port  and  a  gas  tight  gland 
through  which  a  water  cooled  RF  coil  is  passed.  A  high  capacity  diffusion 
vacuum  pump  is  also  attached  through  a  series  of  valves  and  a  liquid  nitrogen 
cold  trap  to  the  stainless  steel  Lank.  A  5  kW  RF  generator  is  connected 
to  the  RF  coil  and  provides  the  heat  to  boil  off  the  copper.  Boiling  off 
the  copper  in  the  presence  of  a  vacuum  is  necessary  both  to  prevent  tungsten 
oxidation  and  to  achieve  copper  boiling  at  relatively  low  temperatures 
(copper  boils  at  2200°F  at  a  pressure  of  2.5  x  10-3  mm  of  Hg) .  This  facility 
was  capable  of  achieving  both  these  conditions. 

B.4  Schlieren  System 

The  Schlieren  system  employed  for  this  study  (Figure  B-4)  was  manu¬ 
factured  by  Baird  Atomic  Inc.  It  consists  of  a  xenon  light  source  that  is 

collected  into  a  parallel  beam  and  passed  over  the  CONAP  test  specimen. 

The  parallel  beam  is  then  focused  and  used  to  develop  an  image  on  either  a 
viewing  screen  or  camera  film.  In  the  portion  of  the  test  section  where 
the  density  is  varying  due  to  helium  injection  the  parallel  light  rays  will 
be  bent  and  consequent. !>  not  pass  through  a  focal  point.  A  knife  edge 

is  inserted  at  the  focal  point  to  eliminate  rays  that  have  been  deflected 

from  parallel.  The  elimination  of  these  rays  from  the  image  results  in  a 
variation  of  illumination  which  is  proportional  to  the  density  gradients 
across  the  helium. 
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Figure  B-3.  Vacuum  Evaporation  Equipment  for 
Removal  of  Copper  from  Tungsten 


Figure  B-4.  Schlieren  System 
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HEAT  TRANSFER  BETWEEN  A  POROUS  MEDIUM,  WHICH  IS 
UNDER  A  UNIFORM  HEATING  MODE, 

AND  A  COOLING  FLUID 
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appendix  c 


Heat  transfer  between  an  elemental 
may  be  written  as  in  (1) 


volume  in  the  porous  tungsten 


q  = 


k 


m 


dT 

m  +  Q'x. 
dx 


(1) 


where 


q  «  heat  transfer  between  elemental  volumes 
km  =  thermal  conductivity  of  matrix 
Tm  =  matrix  temperature 
q'  =  internal  heat  generation 


Differentiating  (1) 
ductivity  ui 
throughout  the  porous  matrix  yields 


with  respect  to  x  and  assuming  that  the  thermocon 


(k  )  as  well  as  the  internal  heat  generated  are  constant 
v  nr 


dT  A  nl 

dq  =  k  m  +  Q  . 
dx  dx2 


(2) 


Heat  may  also  be  exchanged  between  the  porous  tungsten  and  the  transpira¬ 
tion  coolant  by  convection  as  shown  in  Equation  (3) 


where 


h'  -  heat  transfer  coefficient  (volumetric) 
Tc  =  coolant  temperature 


Differentiating  with  respect  to  x  and  assuming  that  the  internal  heat 

transfer  coefficient  h'  is  constant  with  respect  to  x  and  that  h 

>  h’x  d  (T  -T  )  yields 
rr~  m  c 
dx 


dq  =  h '  (T  -T  )  . 

me 


(4) 


Since  the  coolant  temperature  within  the  porous  tungsten  ^expected  t to  1 be 
above  the  critical  temperature  for  NH3  (271°F)  and  N2  (-233  F) ,  the  he 
of  vaporization  was  neglected  and  all  heat  absorbed  as  shown  in 


q  =  m  C  (T  -T.) . 
n  pci 

c 

where 

Cp  -  coolant  specific  heat 
Tjl  «  initial  temperature 
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Assuming  that  the  coolant  specific  heat  remains  constant  (5)  may  be 
differentiated  to  yield 

dq  =  m  C  ^C.  (6) 

pc  dx 

Setting  (4)  equal  to  (6)  and  differentiating  results  in 

h'  -  h'  <3Tc  =  m  C 

dx  dx  Pc  dx^- . 


Differentiating  again  results  in  (7) 


d  T, 


m 


h* 


dx2 


d2Tc 

&T 


m  C 


d3Tc 

dP“ 


(7) 


To  solve  for  either  the  coolant  or  matrix  temperature,  one  or  the 
other  must  be  eliminated  from  (7) .  Since  in  the  steady  state  all  heat 
introducted  into  an  element  must  be  removed  by  the  coolant,  equation  (2) 
may  be  set  equal  to  (6) 

2 

km  d  Tm  +q'  =  m  Cpc  ^Tc . 
dx2  dx 


dl 


c  +  h' 


dx- 


mC 


ir  to  (7 

and  grouping 

r—  -1 

d2Tc 

-  hj_ 

dTc 

dx2 

dx 

=  0. 


Redefining  coefficients  in  (8)  as  shown  below 


C 


(8) 


The  solution  of  Equation  (8)  with  the  simplified  coefficients  shown  above 
is 


Tc  =  C,  +  C2  exp  X{_-A+ 

V 


/-A+Va2  -  4b\  „  x/-A-Va2-4b\  .  „  .. 

\ — 2 — )  +c3exp  V — 2 — ;  +  §x- 


(9) 


To  solve  (9),  the  three  constants  Clf  C2,  and  C3  must  be  evaluated. 

The  evaluation  of  these  three  boundary  conditions  are  based  on  the  following 
three  assumptions!  q' 


1) 


At  the  porous  matrix  backside 
the  coolant  temperature  is  equal  to 
the  coolant  supply  temperature  or 
at  X=0,  Tc=Ti . 


m 


[4-  l  _ 


2) 


At  the  porous  matrix  backside  the  heat  conducted  in  the  coolant  is 
equal  to  the  heat  conducted  ir  the  matrix  or 


at  x  =  0;  k 


m 


dT 
_ 1 

dx 


m  =  k 


dT 


dx 
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3)  The  third  constant  may  be  evaluated  based  on  the  assumption  that 
all  heat  generated  by  the  RF  generator  in  the  porous  matrix  is 
reflected  by  the  increase  in  coolant  temperature 


at  x 


L;  T 


Q'L  +  T.  . 
mCpc 

TL_  above  three  boundary  conditions  result  in  three  equations  which  must 

be  solved  simultaneously.  The  numerical  results  are  then  substituted  into 

Equation  (9)  and  the  coolant  temperature  Tc,  computed  in  terms  of  h',  L, 

m,  Cpc,  k  »  k  ,  T.,  and  Q' . 
m  c  i 


Next  the  matrix  temperature  Tm  may  be  computed  as  a  function  of  Tc 
and  x.  Setting  (4)  equal  to  (6)  and  solving  for  Tm 
C 

Tm 

h'  dx 


’P  dT  , 
m  rc  c '  H-  T 


(10) 


dT 


Solving  for  c  by  differentiating  (9)  yields 

dx  .  _ _  .  .... 

dTC  =  C 2  I  -A+  Va2-4B  \  ^  X  / -A+Va2-^ 


dx 


^-A+  VA2-4B  j  exp  x  ^~A+ V'a2~4B^ 


+C- 


^-A- VA2-4fij  exp  v  ^-A-  Va2-4B^ 


C 

B 


(ID 


Matrix  temperature  Tm  may  now  be  numerically  determined  by  substituting 
(11)  into  U0)  and  solving. 


The  pressure  distribution  in  the  coolant  is  based  on  the  modified 
Darcy 1 s  law 


my  +  gp  /m\  2 

_rP  \f>)  j 


dp  =  -  my  +  gp 
dx  Tp 


(12) 


Where  T  is  the  permeability  and  g  is  the  viscous  resistance  coefficient. 
Substituting  the  gaseous  equation  of  state  for  p  =  P_  into  (12) 
and  integrating  yields  BRT 

n  h 


2  2  2 
P  -  P  =  7  *  RLZ 
o  I  m 


I 


y  Tdq  +  2gm  RLZ  J  Tdn 
o  ° 


(13) 


where  q  =  x. 

L 

The  gas  viscosity  M  is  a  variable  as  a  function  of  local  coolant 
temperature.  This  variation  may  be  approximated  by  (14) 


y. 

y  * 


ft-) 


N 


(14) 


04 


. ‘ . . . .  ■  - - - - - - - 


Where  p*  and  Tc*  and  the  exponential  N  are  constants.  These  constants  are 
shown  below  for  N2  and  NH3 . 


Coolant 
T*  (°R) 

p* (lb/ft  sec) 
N 


Nitrogen 

560 

1.13  x  10-5 
.70 


Ammonia 

540 

.979  x  10-5 
.52 


Substituting  (14)  into  (13) 
Note  that  the  integrals 
in  this  analysis  employing 


f 

Jo 


,  the  local  coolant  pressure  may  be  determined. 


r  n 

T  dr)  and  |  TM+ 
^  o  C 

finite  difference 


dn  were  solved  numerically 
techniques. 
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00010 
00020 
00030 
00031 
00040 
00050 
000G0 
00070 
00080 
00000 
00100 
00110 
00120 
00130 
00140 
00150 
001G0 
00170 
00180 
00100 
00200 
00210 
00220 
00230 
00240 
00250 
00200 
00270 
00280 
00290 
00300 
00310 
00320 
00330 
00340 
00350 
003G0 
00370 
00380 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
00460 
00470 
00480 
004  90 
00500 
00510 
00520 
00530 
00540 


DOUBLE  PRECISION  XIII, Y, CO 
REAL  M  N  LAM 

DIMENSION  P(20),P0(20),M(20),Cl(20),C2(20),XIN(20),Y(20),CO(ll) 
DIMENSION  A1(20),TC1(20) 

1  WR I TE ( 6 , G  0  00 ) 

6000  FORMAT (111  ,47HCOLD  FLOW  CHARACTERISTICS  OF  SINTERED  MATERIALS/) 
WRITE(6,  GOOD 

G001  FORMAT  (111  ,42!IINPUT  0  IF  PRESSURE  DROP  IS  TO  BE  COMPUTEP/lil  ,5GNI! 
1PUT  1  IF  DISK  PERMEABILITY  AND  BATA  ARE  TO  BE  COMPUTED/) 

READ( 5, * )  A 
UR  I TE(6, 6002  ) 

G002  FORMAT ( 111  ) 

I F ( A) 100, 100,200 

100  UR  I TE ( G, GOO  3 ) 

6003  FORMAT  ( 111  ,  D1HINPUT  0  IF  COOLANT  IS  LIQUID  OR  1  IF  COOLANT  IS  GAS/ 
1) 

READ( 5, * )  B 
UR  I TE(C, 6004 ) 

G004  FORMAT  (111  ) 

I F( B ) 10 1, 101, 150 

101  UR  I T E ( 6 ,  6005 ) 

GO  05  FORMAT  ( 111  ,  59UINPUT:  PERMEABILITY  FT2  ,  DENSITY  LP/FT3  , 

1  LB/SEC/1M  , 45HFLUI D  VISCOSITY  LB/FT-SEC  ,  DISK  THICKNESS  IN/lll  ,4 
2GHINERTIAL  RESISTANCE  COEFF  1/FT  ,FLOU  AREA  IN  2/) 

READ( 5, * )  LAM, RHOL, V I SCL,  X, BATA, A 
A-A/144. 

X=X/ 12 . 

700  UR  I  T E  ( 6 G  0 4 0  ) 

G040  FORMAT  (111  ,2 311  INPUT:  MASS  FLOW  LB/SEC/) 

READ ( 5, * )  W 
W-W/A 

DP=((U*VISCL*X)/(LAM*RIIOL*32.174))+(((W**2. )*CATA*X)/(RHOL*32. 174) 
1) 

GO  TO  500 
150  UR  I TE( b,  6006  ) 

GOOG  FORMAT  ( 111  ,  58HHJPUT:  PERMEABILITY  FT2  ,  LR/SEC  ,COMPRESSI 

13  I  L  I  TY/lll  ,  4  2I1GAS  CONSTANT  FT-LB/LB  R  ,GAS  TEMPERATURE  R/lll  ,49111)1 
2SK  THICKNESS  IN  ,  INERTIAL  RESISTANCE  COEFF  1/FT/lll  ,  G 1 H  R  F  F  TEMP  P 
3, REF  VISC  LB/FT-SEC  , V 1  SC  EXPONENT  ,FLOU  AREA  liJ  2/) 

READ (5,*)  LAM, Z, R, TC,  X,  BATA, TR, V I  SCR, N, A 
A-A/144. 

X  =  X / 12  . 

705  WR I TE (6, 6041 ) 

G041  FORMAT ( 111  ,23IIINPUT:  MASS  FLOW  LB/SEC/) 

READ ( 5, * )  W 
U-U/A 

A2  =  ( ( 2 . *W*R*X*Z*V I  SCR) /(  LAI ;*32 . 174 ) ) * ( ( TC** ( N  +  l. ) ) /( ( TR**N) ) ) 

B 1- ( ( 2 . *3ATA*R*X*Z*TC*(W**2 . ) ) /32 . 174 ) 

OP2-A2  +  151 
GO  TO  550 
200  UR  I TE(G,  G007) 

COO  7  FORMAT  (HI  ,  51111 1'PUT  0  IF  COOLANT  IS  LIQUID  OR  1  IF  COOLANT  IS  GAS/ 
1) 

READ( 5, *  )  C 
UR  I TE( G,  6008 ) 


00550 

00560 

00570 

00580 

00590 

00600 

00620 

00630 

00640 

00650 

00660 

00670 

00680 

00690 

00700 

00701 

00710 

00711 

00720 

00730 

00740 

00741 

00750 

00760 

00770 

00780 

00790 

00800 

00810 

00820 

00830 

00840 

00850 

00860 

00870 

00880 

00890 

00900 

00910 

00920 

00930 

00940 

00950 

00960 

00970 

00980 

00990 

01010 

01020 

01030 

01040 

01050 

01060 

01070 

01080 

01090 

01091 

0110- 


6008  FORMAT ( 1H  ) 

IF(B)201, 201,250 

201  WRITE(6,  6009 ) 

6009  FORMAT ( 1H  ,56HINPUT:  LIQUID  DENSITY  LB/FT3  ,  LIQUID  VISCOSITY  LB/FT 
1-SEC/1H  ,  2 1HSPEC I  MEN  THICKNESS  IN/) 

READ( 5, * )  RHOL, V  I SCL, X 
X-X/12. 

205  WRI  TE(6, 6010 ) 

6010  FORMAT ( 1H  ,22HINPUT:  NO  OF  DATA  SETS/) 

READ(5, * )  NNU 

WRITE(6, 6011) 

6011  FORMAT (1H  ) 

WR I TE(6, 6  0 12 ) 

6012  FORMAT ( 1H  ,55IIINPUT:  UPSTREAM  PRESSURE  PSIA  ,  DOWNSTREAM  PRESSURE  P 

IS  I  A/  IN  ,  1311  FLOW  AREA  IM2/1II  ,47HCOOLANT  FLOW  RATE  LB/SEC  (M  IN  ASC 
2ENDIIJG  ORDER/)  / 

READC  5, * )  (P( J ), J  =  l,  NNU), ( PO(J ),  J  =  l, NNU), (Al( J ) ,  J  =  l, NNU), (M( J ) , J= 

11, NNU) 

WRITE(G, 6013) 

6013  FORMAT ( 111  ) 

DO  210  J  =  1 ,  N  N  U 
A1(J)=A1( J)/144. 

P( J )=P( J ) *144 , 

P0( J)*PO( J )*144 . 

M(J)=M(J)/A1(J) 

Cl(J)=((P(J)-POCJ))*(32.179*RHOL))/(VISCL*X*M(J)) 

C2 ( J ) =M( J ) / V  I  SC L 
210  CONTINUE 
NDATA-HNU 
MD-1. 

DO  215  J=1/NNU 
Xlll{  J)  =  C2  ( J  ) 

Y(J)=C1(J) 

215  CONTINUE 

CALL  POLLSQ  (IJDATA,  X  I N, Y, ND, RMS, CO ) 

URITE(6,6014) RMS 

6014  FORMAT  ( 111  ,  4IIRMS  =  1X/  E 10 . 3/ ) 

WRITE(6/6015)CO(1),CO(2) 

6015  FORMAT  (Hi  , 6HC0( 1 ) =,  IX, E 10 . 3/ 1H  , 6MCO( 2  )  =, IX, E 10 . 3/ ) 

GO  TO  600 

250  WRITE(6,G016> 

6016  FORMAT ( 1H  ,  4611 1 NPL'T .  COMPRESS  I  B I  L I TY  ,GAS  CONSTANT  I"T-LB/LB-R/1H  ,3 

16H  /DISK  THICKNESS  IN/lll  ,  66HREFERENCE  TEMP  P  .RE 

2F  VISC  LB/FT-SEC  ,VISC  EXPONENT  ,  /) 

READ( 5,  * )  Z,fi,X,TR,VISCR,N 
WR  I TE ( 6, GO  19 ) 

6019  FORMAT (  Hi  ) 

X=X/12 . 

255  WRITECG,  G017) 

6017  FORMAT  ( 111  ,22HINPUT:  NO  OF  DATA  SETS/) 

READ( 5,  * )  NNU 

WR I TE( 6 ,  G  0 18  ) 

6018  FORMAT  (111  ) 

WRITECG,  6020) 

6020  FORMAT ( 111  ,55HINPUT:  UPSTREAM  PRESSURE  PSIA  , DOUNSTPEAM  PRESSUP.E  P 
1SIA/1II  ,  32HGAS  TEMPERATURE  R  ,FLOU  AREA  IN2/1H  .  4  7HCOOLANT  FLOW  RA 
2TE  LB/SEC  (M  IN  ASCENDING  ORDER)/) 

REAIK 5,  *  )  (  P(  J ) ,  J  =  1,  NNU) /  ( PO(  J) ,  0  =  1,  IJNU) ,  (TC1(  J  ) ,  0  =  1,  NNU) ,  (Al(  J  ) , 
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PI pfjj? 


01101 

1J«1,NNU),  (IKJ),  J=1,NNU) 

OHIO 

WRITE(6, 6021) 

01120 

6021 

FORT  '.AT  ( 1H  ) 

01130 

DO  260  J  =  1,IINU 

01131 

Al( J ) =A1 ( J )  / 144 . 

01140 

P(J)-P( J  )  *  144 . 

01150 

P0(  J  )  =  PO  (  J  )  *  144  . 

01160 

M(J)-M(J)/A1(J) 

01170 

Cl( J)«( (  32.174*(P(J)**2.-PO(J)**2.  )))/(( 2 . *V 1 SCP*M(J ) *Z*R*TC1  ( J ) *X 

01130 

l*((TCl(J)/TR)**tJ))) 

01190 

C2(J)=M(J)/(VISCR/((TC1(J)/TR)**H)) 

01200 

260 

COUTirJUE 

01210 

U  DAT  A  =  f  J  M  U 

01220 

HD-1. 

01230 

DO  265  J-l, NNU 

01240 

X 1 M ( J)»C2(J) 

01250 

Y(J)=C1(J) 

01200 

265 

CONTINUE 

01270 

CALL  POLLSQ  ( NCATA, X 1 N, Y, HD, RMS,  CO ) 

01280 

UR  1 TE ( 6, 0022 ) RMS 

01290 

6022 

FORMAT ( 111  ,  4MRUS-1X ,  E 10 . 3/ } 

01300 

URITE<G,G023)CO(1),CO(2) 

01310 

6023 

FORI’AT  ( 111  ,GHCO(l)=,  IX,  E10. 3/111  ,  6HCO(  2  )  =  ,  IX,  E 10 . 3/ ) 

01320 

60  TO  600 

01330 

500 

UR  1 TE ( G , G024 ) DP 

01340 

6024 

FOR? '.AT (11!  ,  2 1!  II.  1 QU 1  D  PRESSURE  DROP-, IX, E10 . 3,  6HLR/FT2/ ) 

01341 

UR  1 TE (6,  6045 ) 

01342 

6045 

FORMAT  ( 1H  ,  G4HI IJPUT  0  IF  PROBLEf;  IS  TO  BE  CONTINUED  OR  1  IF  PROCLE 

01343 

1M  IS  CHANGED/) 

01344 

READ(  5,  *  )  I; 

01345 

UR  1 TE ( 6,  6  046  ) 

013  4  G 

G  046 

format  (111  ) 

01347 

1 F ( n  >  700,  700,  1000 

01300 

5  50 

UR  1 TE ( 6 , G  02  5  ) DP2 

01370 

6025 

FORMAT ( 111  , 2GHGAS  PRLSSURC  DROP  SQUARED-, IX, F10. 3, 7HLB2/FT4/) 

01371 

WRITE (6, 6047) 

•j  13  72 

6047 

FORMAT  ( Hi  ,64111 NPUT  0  IF  PROCLEM  IS  TO  BE  CONTINUED  OR  1  IF  PR03LE 

01373 

111  IS  CHANGED/) 

01374 

REALM  5,*)  9 

01375 

UR  1 TE (G, C  J-8 ) 

01370 

6043 

FORMAT ( 111  j 

01377 

1  F( D)  705,  705,  1900 

01390 

600 

L  All-1. /CO(l) 

01400 

J-NNU 

01410 

Y2=CO( l)+(CO(2 )*C2 ( J ) ) 

01420 

I5ATA- ( Y2-CO(  1)  )/C2  (J) 

01430 

DO  605  J  - 1,  NNIJ 

01440 

UR  1  T  E  (  6  ,  6  0  2  6  )  C  2  ( J  )  ,  C 1  ( J  ) 

01450 

6026 

FORMAT  (111  ,  71  IP/ VI  SC-,  IX,  E 10. 3/ 111  ,31II(G)(DP)/2(VISC)(M)(Z)(P)(T)(X 

0140  0 

1)=, IX, E1Q . 3/ ) 

0  14  79 

605 

CONTINUE 

OI400 

UR  ITE (6,  6027  ) LAM,  BATA 

01490 

6027 

FORMAT  ( 111  ,1311  PER* '.FACILITY-,  IX, CIO. 3, 3IIFT2/1H  ,  37H 1 MEPT 1  M.  RESISTA 

01500 

INCH  RESISTANCE  COE  :F  =  ,  IX, E10, 3, 4H1/FT/ ) 

01510 

1000 

GO  TO  1 

01520 

1001 

STOP 

01530 

END 
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APPENDIX  E 


A  FORTRAN  IV  LISTING  FOR  THE 
DETERMINATION  OF  THE 
INTERNAL  HEAT  TRANSFER  COEFFICIENT 


E-1 


00010  REAL  N,MDOT,KM,KF,L,MU,  I 

00020  DIMENSION  ANU ( 20 ) , TF ( 20 ) , TM( 20 ) 

00030  1  UR  I TE (6,  GOOO ) 

00040  6000  FORMAT  ( i:  i  .65IIONE  D I  MENS  I  ONAL  TRANSPIRATION  COOLING  MODEL  WITH  I  NT 

00050  1ERNAL  HEAT  HG/) 

00060  HR  I TE (6, 6  jO  3 ) 

00070  6003  FORMAT  ( 111  ,  6911 1 NPIJT:  K  MATRIX  BTIJ/FT  SEC  R,K  FLUID  BTU/FT  SEC  R,CP 

00080  1COOLANT  I3TU/LB  R/) 

00000  READ( 5, *)  KM,KF,CP 

00100  HR  I TE(6, G004  ) 

00110  6004  FORMAT ( 111  ) 

00120  HR ITE(G,6005) 

00130  6005  FORMAT  ( 111  ,G0!HNPUT:Q  INTERNAL  HEATING  BTU  /FT-CU  SEC, COOLANT  INLE 

00140  IT  TEMP  R/ ) 

00150  READ ( 5, * )  Q I N, TF  I 

00160  HR ! TEC  6, G  006 ) 

00170  600G  FORMAT ( 111  ) 

00180  UR  I TE ( G, 6 00 7  ) 

00190  6007  FORMAT  ( 1H  ,  5411 1  NPIJT:  GAS  CONSTANT  FT-LB/LB  R,  COMPRESS  I  B  I  L  I  TV  IJIIITLE 

00200  1SS/1U  , 23HLIQUI D  DENSITY  LE/FT  CU/) 

002  10  READC  5,  * )  R,Z,RIIOL 

00220  UR  I TECG, G008  ) 

00230  6008  FORMAT ( HI  ) 

00240  UR  I TE (6,  6  009 ) 

002  50  6009  FORMAT  ( 111  ,  65111  NPIJT:  FLUI  D  VISCOSITY  LB/FT  SEC, GAS  PCFERENCE  VISCOS 

00260  1ITY  LB/FT  SEC/lil  ,29HGAS  REFERENCE  TEMP.  R,VIS  EXP/) 

00270  READC  5, * )  MU,  V  I  SC, TEMP, N 

00280  UR  I TEC6, G010 ) 

00290  6010  FORMAT ( l!l  ) 

00300  40  UR  I TECG, 6001 ) 

00310  6001  FORMAT  ( 111  ,  6  5HI f'  PUT :  MOOT  LB/FT  SO  SEC,  INTERNAL  HEAT  TRANS  COEFF  CT 

00320  1U/FT  CU  SEC  R/lll  ,5411  HA  L  L  THICKNESS  I  NCMES,  MATR I  X  PERMEABILITY  F 

00330  2T  SQ,  TU  R/lll  ,  3011 1  NERT  I  AL  RESISTANCE  COEFF  1/FT/) 

00340  READC  5,  * )  MDOT,  IIPR,  L,  GAM,  TU,  BATA 

00350  UR ITE(6,6025) 

00360  6025  FORMATC1II  ) 

00370  UR  I TECG, 6011) 

00380  6011  FOP.MATC  111  ,  4  511 1  NPUT :  NO  OF  HALL  SEGMENTS,  NORMAL  I  7.  CD  LOCATION/) 

00390  READC  5, * )  NNU, ( ANUC J ) , J  =  l, NNU) 

00400  UR  I TE ( 6,  6 0 12 ) 

00410  6012  FORMAT ( 111  ) 

00420  C 

00430  L  =  L / 1 2 . 

00440  A  =  HPR/(MDOT  *CP) 

00450  C  =  -  (IIPR /KM) 

00460  70  C  =  (IIPR*QI  II )  /  C MDOT* CP* KM) 

00470  D=(KM*MDOT*CP)/HPR 

00480  li=(-A  +  SERT(A**2.-(4.*B)  ))/2. 

00490  l=(-A-SQRT(A**2.-(4.*R)))/2. 

00500  E  =  ( D* I )  +  KM  +  KF 

00510  F=KM-KF 

00520  G  =  KF  —  KM— ( H  *  D ) 

0  05  30  ANUM=  ( ( C*F )  /  (  B*II*G ))-(((  C*F )/( B*II*G ) ) *EXP(  L*H ) )  +  ( C 0  I  N* L)  /  (MDOT* CP ) 

00540  l)-( ( C*  L ) /D ) 

0  0550  ADNOM=  ( ( (  I  *E  )  /  ( G*H ) )  *EXP(  L*ll ) )  + ( EXPC  L*  I  )  )  -  ( (  I  *F  )  /  ( G*H )  )  - 1 . 

00560  C3=ANUM/ADNOM 

0  05  70  C  2  =  ( C  C  3  *  I  *E)/(H*G) )+( (C*F)/(B*H*G) ) 

00580  C1=TF I -C3-C2 

°°590  URITE(6,  6020)L,A,B,0,D,!I,  I ,  E,  F,  G,  ANIJM,  ADNOM,  C 1,  C2  ,  C3 


E-2 


c 


OOGOG 
OOGIO 
00G20 
00630 
00640 
00650 
00660 
00670 
00680  C 
00600  C 
00700  C 
007  10 
00720 
00721 
00730 
00740 
00750 
00760 
00761 
00762 
00763 
00764 
00765 
00766 
00767 
00768 
00769 
00770 
00730 
00810 
00811 
00820 
00830 
00840 
00850 
00860 
00870 
00871 
00872 
00880 
00881 
00382 
00890 
6  09  00 
00910 
00920 
0  03  30 
00940 
00950 
00360 
00970 
00380 
00990 
01000 
OlOlO 
01020 
01030 
01040 
01050 
01060 


6020  FORMATUIl  ,  2IIL-E10 . 3/ 111  ,2HA-E10.3/lh  2HP.»,  E 10 . 3/ 1H  /  2IIC  =  ,  Fl°  . 3/ 1 

111  ,2111)=, CIO.  3/111  ,  2MII  =  ,  E10. 3/111  ,2111  =,tlO.  3/11!  ,  2,l£"'  E  3'  2 1, 

2=  E10. 3/Hi  ,  2 i IG  = ,  E10 . 3/111  ,  5HANUM*,  E10 . 3/lil  , Cl lAPMOf E  10  .  3 / 111  ,3 
3C1  =  ,  E10. 3/111  ,  3HC2-,  E10. 3/1H  ,  3IIC3-,  E10 . 3/ ) 


URITE(G,60I3) 
6013  FORMAT ( 2 X  ,65KX 
1  I) PIJ2 G /  ) 


If,  R 


TF  R 


DPN2 


CHECK  FOR  GAS  OR  FLUIP 


I  F  (  RIIOL .  GT .  0 .  )G0  TO  500 

5uri=o.  o 
suni=o.o 

DO  200  J  =  1, HMU 
JJ  =  J 

XaL*ANU( J )  ,,„v 

TF(,J)  =  (C2*EXP(X*II))  +  (C3*CXP(X*I  )  )  +  ((C*X)/B)*Cl 
DD=TF(J)-TFI 
I  F  (  619)100,  100,  101 

100  TF ( J ) =TF I 

BDT  =  (  (C2*II*EXP(X*II)  ) )+ (C/B ) 

TM( 0 ) = ( 6DT/A) +TF ( J ) 

CC  =  TM( J ) -TF  I 
I F ( CC)  105,  105,  102 
105  TI1( J )  =TF I 

101  DTF=((C2*II*EXP(X*II)))+(C3*I*EXP(X*I  ))  +  (C/3) 

TU(  J)-0)TF/A)*TF(  J) 

102  COIJS=333580. 032 
CONS1=1334320. 1 


dp;j2=((r*hgot*i;u*l*z)/(coijs*ga;‘))*c,ji- 
I F ( J . EQ. 1)DPN2S=0 . 

30  |l|NT=>('(l(TF(  JjI**(‘1  +  1)  +  TF(  JJ--l)**(ri+l)))/2.  )*(ANU(  JJ)-ANU(tlJ-l) ) 

T VfIT  1=  ( < TF ( JJ)+TF(  JJ-1) )  / 2  .  )  *  ( AIJIK  J  J)-AtJU(»l  J-l ) ) 

suni=suni+Ti:iTi  ,  MXWCltf, 

n P f J 2 3 1  =  ( 1 .  /  ( C0!i5*GAM) )  *'T0  T*R* L*7.* V I  SC* ( TEMP**  ( -N) )  *SUf 
!)PN252=(2./CONSl)*nATA*(r:nOT**2.  ) *R*L*7.*SUt!l 
DPn2S  =  DPII2Sl  +  DPN2S2 
20  COUTINUE 

UR^TL(6/6014)X,TII(J),TF(J),DPf!2/DPM2S 

6014  FORMAT ( 5E14 . 5 ) 

200  CONTINUE 
GO  TO  910 


C 


C  FLUID 

C 

500  CODT I  HUE 

UR  I  TC(G,G021)  . 

60  21  FORMAT  (111  ,  2  7IILIQUIR  OPTION  NOT  AVAILABLE/) 

GO  TO  1000 

900  do  600  J  =  i,unu 

X  =  L*AI!IJ(J) 

f>E!!OM=E*Rl*EXP(Rl)  +  F*R2*EXP(R2) 

BP.KT  =  (  C*EXP(  R1*ANU(  J  )  )+F*EXP(R2* ANU( J ) ) )  /DEN0M 
Tll(  J)=B*C*CRKT+TF1 
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J1070 

RKT*  ( E*R1*EXP(R 1*ANU( J ) )+F*P2*EXP(R2*ANU( J ) ) ) /DENOM 

01080 

TF( J)=C*RKT+TF 1 

01090 

DPH*( ( AMDT*AMU*L) / (4G33. 056*GAM*RH0L) )*ANU( J ) 

01100 

X=X*12 . 

OHIO 

URITE(6,6015)X,TM/TF(J),DPfJ 

01120 

6015 

FORMAT ( 5E15 . 6 ) 

01130 

600 

CONTINUE 

01140 

1000 

CONTINUE 

01150 

RR= ( MDOT  *CP) / L 

01100 

NIJ-1 

01161 

DD= .  5 

01170 

50 

DT=DD*(TW-TM( NN) ) 

01180 

1 FCDT.GT.5.  )GO  TO  700 

01190 

DT»DD*(TM(NN)-TW) 

01200 

1  F ( l)T.  GT.  5 . ) GO  TO  800 

01210 

60 

WR 1 TE ( 6, G  017 ) Q 1 M 

01220 

6017 

FORMAT ( 111  ,  17H 1  MTERNAL  HEAT  1  NG«,  IX,  E10 . 3,  13IIBTU/FT 

01230 

GO  TO  40 

01240 

700 

Q 1  f!  =  RR*  ( (TF  (  UN  )  -  TF  1  )  +  (DT  > ) 

01250 

GO  TO  70 

012CC 

800 

Q  1  H  =  RR* ( ( TF  (  NN ) -TF I )-(DT)) 

01270 

GO  TO  70 

01280 

910 

GO  TO  1 

01290 

90 

STOP 

01300 

END 
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